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We study the abundance of a particle species in a thermalized plasma by introducing a quantum 
kinetic description based on the non-equilibrium effective action. A stochastic interpretation of 
quantum kinetics in terms of a Langevin equation emerges naturally. We consider a particle species 
that is stable in the vacuum and interacts with heavier particles that constitute a thermal bath in 
equilibrium. Asymptotic theory suggests a definition of a fully renormalized single particle distri- 
bution function. Its real time dynamics is completely determined by the non-equilibrium effective 
action which furnishes a Dyson-like resummation of the perturbative expansion. The distribution 
function reaches thermal equilibrium on a time scale ~ l/2rfe(r) with rfe(r) being the quasiparti- 
cle relaxation rate. The equilibrium distribution function depends on the full spectral density as a 
consequence the fluctuation-dissipation relation. Such dependence leads to off-shell contributions to 
the particle abundance. A specific model of a bosonic field $ in interaction with two heavier bosonic 
fields xi,2 is studied. The decay of the heaviest particle and its recombination lead to a width of the 
spectral function for the particle "1> and to off-shell corrections to the abundance. We find substantial 
departures from the Bose-Einstein result both in the high temperature and the low temperature but 
high momentum region. In the latter the abundance is exponentially suppressed but larger than 
the Bose-Einstein result. We obtain the Boltzmann equation in renormalized perturbation theory 
and highlight the origin of the differences. Cosmological consequences are discussed: we argue that 
the corrections to the abundance of cold dark matter candidates are observationally negligible and 
that recombination erases any possible spectral distortions of the CMB. However we expect that 
the enhancement at high temperature may be important for baryogenesis. 



I. INTRODUCTION 



Phenomena out of equilibrium played a fundamental role in the early Universe: during phase transitions, baryo- 
genesis, nucleosynthesis, recombination, particle production, annihilation and freeze out of relic particles, some of 
which could be dark matter candidates [lj,|3,|3|. Of the many different non-equilibrium processes, particle production, 
annihilation and freeze-out and baryogenesis^ are non-equilibrium kinetic processes which are mainly studied via 
the Boltzmann equation^, |3, Q . 

The Boltzmann kinetic equation is also the main approach to study equilibration, thermalization and abundance of 
a species in a plasma. A thorough formulation of semiclassical kinetic theory in an expanding Friedmann-Robertson- 
Walker cosmology is given in ref . . 

However the Boltzmann equation is a classical equation for the distribution function with an inhomogeneity deter- 
mined by collision terms which are computed with the S-matrix formulation of quantum field theory. The collision 
term in the Boltzmann equation is obtained from the transition probability per unit time extracted from the asymptotic 
long time limit of the transition matrix element. This is tantamount to implementing Fermi's golden rule. Potential 
quantum interference and memory effects are completely ignored in this approach. Furthermore a single particle 
distribution function, the main ingredient in the Boltzmann equation, is usually defined via some coarse graining 
procedure. All of these shortcomings of the usual semiclassical Boltzmann equations when extrapolated to the realm 
of temperatures and density in the Early Universe, suggest that in order to provide a reliable understanding of such 
delicate processes such as baryo and leptogenesis a full quantum field theory treatment of kinetics may be required^. 

One of the basic predictions of the Boltzmann equation is that the local thermodynamic equilibrium solution for 
the abundance of a particle species is determined by the Bose-Einstein or Fermi-Dirac distribution functions, hence 
exponentially suppressed at low temperatures (in absence of a chemical potential) . 

This basic prediction has recently been challenged in a series of articles 5] wherein a surprising result is obtained: 
the abundance of heavy particles with masses much larger than the temperature is not exponentially suppressed as 
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the Boltzmann equation predicts but the suppression is a power law. Such result, if correct, can have important 
consequences for the rehc abundance of cold dark matter candidates. 

This result, however, has been criticized and scrutinized in detail by several authors^ 0, 1^ who concluded that 
it is a consequence of the definition of the particle number introduced in rcf. . The definition of the total number 
of particles proposed in 5] is based on the non-interacting Hamiltonian for the heavy particle divided by its mass 
plus countcrtcrms, which purportedly account for renormalization effects. The results of references^, 0, Q point 
out the inherent ambiguity in separating the contribution to the energy density from the particle and that of the 
bath and the interaction. The ambiguity in the separation of the different contributions to the energy has been 
studied thoroughly in these references in particular exactly solvable modelsQ, effective field theory 7] or a consistent 
treatment of renormalization effects]^. 

Understanding the limitations of and corrections to the Boltzmann kinetic description and potential departures 
from the predicted abundances is important for a deeper assessment of possible mechanisms of baryogenesis as well 
as for the relic abundance of cold dark matter candidates. In the case of baryogenesis, the applicability and reliability 
of Boltzmann kinetics in the conditions of temperature and density that prevailed in the early Universe warrants a 
critical reassessment 1]. Refinements of the usual Boltzmann equation have been proposed in the literature'^. 

While the work in refs.^, _7, has clarified the shortcomings of the definition of the total particle number pro- 
posed inQ explaining the origin of the power law suppression as a consequence of the ambiguity in this definition, 
what is missing from this discussion is a suitable definition of a distribution function and its real time evolution. 
The Boltzmann equation is a local differential equation that determines the dynamics of the single particle distribu- 
tion function. Therefore in order to clearly assess potential corrections to the equilibrium solutions of the familiar 
Boltzmann equation a suitable distribution function and its dynamical evolution must be understood. 

The definition of the distribution function both in non-relativistic many body theorvfioj as well as in relativistic 
quantum field theorv [Tll Il2j | is typically based on a Wigner transform of a two point correlation function, which is 
not manifestly positive semidefinite. Usual derivations of the Boltzmann kinetic equation invoke gradient expansions 
or quasiparticle (on-shell) approximations which lead to Markovian dynamics. Alternative derivations of the kinetic 
equations'13^ which explicitly implement real time perturbation theory often invoke a long time limit and Fermi's 
Golden rule which enforces energy conservation in the kinetic equation. This is also the case in the dynamical 
renormalization group approach to quantum kinetics advocated in ref.0| although this latter method allows one to 
systematically include off-shell corrections. Whichever method of derivation of the kinetic equation is used, the first 
step is to define a single particle distribution function. 

Any definition of the distribution function of particles that decay in the vacuum (resonances) is fraught with 
ambiguities because the spectral representation of such particles is not a sharp delta function but typically a Breit- 
Wigner distribution. Since these particles decay even in vacuum and do not exist as asymptotic states any definition 
of an operator that "counts" these particles will unavoidably be ambiguous. 

In this work, we circumvent this ambiguity by focusing on the study of the quantum kinetics and equilibration 
dynamics of the distribution functions of particles that are stable at zero temperature associated with a field 
Stable physical particles are asymptotic states which can be measured and a distribution function for the single 
particle physical states can be introduced according to the basic assumptions of asymptotic theory. While our 
ultimate goal is to find a quantum kinetic description for phenomena in the early Universe, in this article we focus on 
a study in Minkowski space-time as a first step towards that goal. 

Goals and methods: 

In this article we provide a framework for non-equilibrium quantum kinetics beyond the usual Boltzmann equation. 
This non-equilibrium formulation includes off-shell and non-Markovian (memory) processes which are not accounted 
for in the semiclassical Boltzmann equation and result in modifications of the equilibrium abundances. We focus on 
the case of a scalar field $ coupled to other heavier fields for a wide variety of relevant interacting quantum field 
theories. Here we consider that the heavier fields constitute a thermal bath in equilibrium. In order to study the 
thermalization of the $ particle as well as the time evolution of its distribution function we consider the case in which 
the field $ is coupled to the thermal bath at some initial time t^. We then obtain the non- equilibrium effective action 
for the field $ by integrating out the degrees of freedom of the thermal bath to lowest order in the coupling of the 
field to the heavy sector but in principle to all orders in the couplings of the heavy fields amongst themselves. 

At zero temperature the ^-particles are stable because they are the lightest, therefore they are manifest as asymp- 
totic states. Hence according to asymptotic theory we introduce a definition of an interpolating number operator 
that counts these particles, for example as those measured by a detector in a collision experiment in the vacuum. At 
finite temperature the distribution function is the expectation value of this interpolating operator in the statistical 
ensemble. The real time evolution of this distribution function is completely determined by the non-equilibrium 
effective action and its asymptotic long time limit determines the abundance of the physical particles $ in the ther- 
mal plasma. The non-equilibrium approach introduced here, borrows from the seminal work on quantum Brownian 
motion^ Ellllll which is adapted to quantum field theory. 
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After the discussion of the general case, we introduce a specific model in which the scalar field $ associated with 
the stable particle couples to two heavier bosonic fields which constitute the thermal bath. At lowest order in the 
coupling we find that the $ particle despite being the lightest, acquires a width in the medium as a consequence of the 
two body decay of the heavier particle and its recombination in the plasma. These processes result in a broadening 
of its spectral function and corrections to its equilibrium abundance. 

Brief summary of results: 

• We obtain the non-equilibrium effective action for a field <i> coupled to other heavier fields by integrating out 
the latter to lowest order in their coupling to the field but in principle to all orders in the couplings amongst 
themselves. The heavy fields are taken to be in thermal equilibrium and therefore provide a thermal "bath" 
for the $ field. The resulting non-equilibrium effective action can be interpreted as a generating functional of a 
stochastic field theory in which the (integrated out) heavy fields introduce a Gaussian but colored noise and a 
non-Markovian self-energy (dissipative) kernel. 

• We introduce a definition of the single particle distribution function in the general case of a particle that is 
stable in the vacuum. Stable physical particles are asymptotic states which can be measured by a detector. In 
accordance with the results of asymptotic theory, we introduce a fully renormalized interpolating number oper- 
ator whose expectation value in the non-equilibrium state (density matrix) is identified with the single particle 
distribution function. The time evolution of this distribution function is determined by the non-equilibrium 
effective action and is completely specified by the solution of a stochastic Langevin equation with a memory 
kernel and a Gaussian stochastic noise. The properties of the memory kernel are related to the spectrum of 
the noise by a generalized fluctuation dissipation relation. We argue that the time evolution of the distribution 
function is a result of a Dyson resummation of the perturbative expansion provided by the non-equilibrium 
effective action. The single particle distribution function becomes insensitive to the initial conditions at time 
scales longer than the "quasiparticle" relaxation time and its asymptotic long time limit describes a thermalized 
state. 

• A specific example is studied in detail. This is a model of Bosonic scalar fields with a coupling 5 $ xi X2 with the 
masses of the "bath" fields xi,2 obeying the hierarchy Mi > M2 » In this case the particles associated 
with the field $ are stable in the vacuum. However, at finite temperature the particle $ acquires a width from 
the two-body decay and recombination process xi + X2 ■ We study the approach to thermal equilibrium of 
the single $ particle distribution function whose asymptotic long time limit yields their equilibrium abundance 
in the bath. We find that the equilibrium abundance is always larger than that predicted by the Bose-Einstein 
distribution. The enhancement is more significant at high temperatures, as well as at low temperatures but 
large momenta. The departure from the Bose-Einstein result is a distinct consequence of off-shell support of the 
spectral bmction of the $ field in the plasma. 

• We derive the usual quantum kinetic Boltzmann equation in renormalized perturbation theory up to the same 
order in the coupling to the bath as the non-equilibrium effective action. This derivation highlights the neglect 
of memory and correlations in the usual Boltzmann equation. We contrast its prediction for the equilibrium 
abundance, the usual Bose-Einstein distribution, to that from the full quantum kinetic equation with memory 
and off-shell contributions. This direct comparison leads to the conclusion that memory and off-shell phenom- 
ena result in substantial corrections to the equilibrium abundances that are not captured by the Boltzmann 
equation. 

• We conclude that potential corrections to the abundance of cold dark matter candidates as well as distortions of 
the cosmic microwave background post recombination are negligible observationally, but substantial corrections 
in a high temperature plasma may be important for baryogenesis. 

The article is organized as follows: in section (|IH) we introduce the general form of the interacting quantum field 
theories considered and develop the formulation in terms of the non-equilibrium effective action. The effective action 
is obtained to lowest order in the coupling of the field $ to the heavier fields (the bath) and in principle to all 
orders in the coupling of the bath fields amongst themselves. We show that a stochastic formulation in terms of a 
Langevin equation emerges naturally. In section (|III|I we introduce the definition of the fully renormalized interpolating 
number operator and the single particle distribution function based on asymptotic theory. The time evolution of this 
distribution function is completely determined by the solution of the stochastic Langevin equation. 

In section (jIV|l we study a specific model in which the $ field is coupled to two heavy scalar fields with a cou- 
pling 5$Xi X2- This interacting quantum field theory provides an excellent testing ground and highlights the main 
conceptual results. We study the dynamics of the distribution function for the $ particle up to one loop order. The 
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asymptotic distribution function is studied for a wide range of parameters alfowing to extract fairly general conclu- 
sions whose validity goes beyond this specific model. In particular we analyze in detail how off-shell effects result in 
large corrections to the usual Bose-Einstein equilibrium abundance. In section (0 we obtain the usual Boltzmann 
quantum kinetic equation and highlight the main assumptions implicit in its derivation. We contrast the predictions 
for the asymptotic abundance between the non-equilibrium kinetic formulation and that of the usual quantum kinetic 
Boltzmann equation, highlighting that memory and off-shell effects are responsible for the differences in the predic- 
tions. Our conclusions and a discussion on the cosmological consequences are presented in section (|VI|I . An appendix 
is devoted to the explicit calculation of the self-energy in the specific example studied. 

II. GENERAL FORMULATION: THE NON-EQUILIBRIUM EFFECTIVE ACTION 

We focus on the description of the dynamics of the relaxation of the occupation number of a scalar field $ which 
is in interaction with other fields either fermionic or bosonic, collectively written as Xi-i with a Lagrangian density of 
the form 

£[<^>ix),x{^)]=U^mx)]+C^,[x^{x)]+g<^>0[x^{x)] (2.1) 

where 0[xi\ stands for an operator non- linear in the fields Xi and £0.* is the free field Lagrangian density for the 
field <& but C^.[xi{x)] is the full Lagrangian for the fields x including interactions amongst themselves. This general 
form describes several relevant cases: 

• Interacting scalars, for example the linear sigma model in the broken symmetry phase. The interaction between 
the massive scalar and the Goldstone bosons is of the form cttt'^ . In this article we focus on the case of a trilinear 
interaction of the form dijXiXj where the fields xi.2 have masses larger than that of the $ field. 

• A Yukawa theory with x being fermionic fields and $ a scalar field, with interaction This could be 
generalized to a chiral model. 

• A gauge theory in which $ is the gauge field and x is either a complex scalar or fermion fields, the interaction 
being of the form with being a bilinear of the fields. In particular this approach has been recently 
implemented to study photon production from a quark gluon plasma in local thermal equilibrium 1^. This case 
is particularly relevant for assessing potential distortions in the spectrum of the cosmic microwave background. 

• Another possible realization of this situation could be the case in which $ is a neutrino field in interaction with 
leptons and (or) quarks which constitute a thermal or dense plasma. 

• The case of a self-interacting scalar field in which one mode say with wave vector k is singled out as the "system" 
and the other modes are treated as a "bath" . 

In all of these cases the fields Xi a-re treated as a bath in equilibrium assuming that the bath fields are sufficiently 
strongly coupled so as to guarantee their thermal equihbration. These fields will be "integrated out" yielding a reduced 
density matrix for the field $ in terms of an effective real-time functional, known as the inffuence functional^^ in 
the theory of quantum brownian motion. The reduced density matrix can be represented by a path integral in terms 
of the non-equilibrium effective action that includes the inffuence functional. This method has been used extensively 
to study quantum brownian motion llTl llS] and for preliminary studies of quantum kinetics in the simpler case of 
a particle coupled linearly to a bath of harmonic oscillatorSj5i^^20]. 

The models can be generalized further by considering that the interaction between $ and x is also polynomial in 
$. However, in this article we will consider the simpler case described by (|2.1|l since it describes a broad range of 
physically relevant cases, and as will be discussed below this case already reveals a wealth of novel phenomena. As 
we will discuss in detail below most of the relevant phenomena can be highlighted within this wide variety of models 
and most of the results will be seen to be fairly general. 

The relaxation of the distribution function is an initial value problem, therefore we propose the initial density 
matrix at a time ti to be of the form 

p{U) = P<s>,i «) Pxi,i (2.2) 

The initial density matrix of the Xi fields will be taken to describe state in thermal equilibrium at a temperature 
T = 1/13, namely 



(2.3) 
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where H^.{xi) is the Hamiltonian for the fields Xi- We will now refer collectively to the set of fields Xi simply as x 
to avoid cluttering of indices. 

In the field basis the matrix elements of p^^i are given by 

($1/5$,,!$') =p*,.(*;*') (2.4) 

The density matrix for $ will represent an initial out of equilibrium state. 

The physical situation described by this initial state is that of a field (or fields) in thermal equilibrium at a 
temperature T = 1//3, namely a heat bath, which is put in contact with another system, here represented by the field 
Once the system and bath are put in contact their mutual interaction will eventually lead to a state of thermal 
equilibrium. The goal is to study the relaxation of the field <I> towards equilibrium with the "bath" . The initial density 
matrix of the field $ will describe a state with few quanta (or the vacuum) initially. 

The real time evolution of this initial uncorrelated state will introduce transient evolution, however the long time 
behavior will be insensitive to this initial transient. Furthermore, we point out that it is important to study the initial 
transient stage for the following reason. As a particle $ propagates in the medium it will be screened or dressed by 
the excitations in the medium and it will propagate as a "quasiparticle" . Its distribution function will be shown to 
become insensitive to the initial conditions on time scales larger than the "quasiparticle" relaxation time. 

The strategy is to integrate out the x fields therefore obtaining the reduced time dependent density matrix for the 
field $, and the non-equilibrium influence functional for this field. Once we obtain the reduced density matrix for the 
field we can compute expectation values or correlation functions of this field. We will focus on studying the time 
evolution of the distribution function, or particle number to be defined below. 

The time evolution of the initial density matrix is given by 



p{tf) = e~'"^'f-'^^piU)e'"^'f-'^^ (2.5) 
Where the total Hamiltonian H is given by 

i7 = i/$($)+iJ^(x)+i//(<i>,x) (2.6) 

The calculation of correlation functions is facilitated by introducing currents coupled to the different fields. Fur- 
thermore since each time evolution operator in cqn. 12.5|l will be represented as a path integral, we introduce different 
sources for forward and backward time evolution operators, referred to as J"*", J~ respectively. The forward and 
backward time evolution operators in presence of sources are U{tf,ti] J+), U~^{tf,ti, J~) respectively. 

We will only study correlation functions of the $ field, therefore we carry out the trace over the x degrees of 
freedom. Since the currents allow us to obtain the correlation functions for any arbitrary time by simple variational 
derivatives with respect to these sources, we can take — > cxd without loss of generality. 

The non-equilibrium generating functional is given by 



Z[] + ,3-] = TvU{^,U;J+)p{U)U-\^,U,J~) (2.7) 

Where J* stand collectively for all the sources coupled to different fields. Functional derivatives with respect to the 
sources J+ generate the time ordered correlation functions, those with respect to J~ generate the anti-time ordered 
correlation functions and mixed functional derivatives with respect to J"*", J~ generate mixed correlation functions. 
Each one of the time evolution operators in the generating functional (|2.7() can be written in terms of a path integral: 
the time evolution operator U{oo,ti; J+) involves a path integral forward in time from to t = oo in presence of 
sources , while the inverse time evolution operator U~^{oo, ti, J~) involves a path integral backwards in time from 
t = oo back to ti in presence of sources J~ . Finally the equilibrium density matrix for the bath e~^'^^ can be written 
as a path integral along imaginary time with sources . Therefore the path integral form of the generating functional 
(1^ is given by 



with the boundary conditions $+(x, i^) = ^i{x) ; <f> {x,ti) — $j(a;). 
The non-equilibrium action is given by 



/oo 
dtd^x [£o,*($+) + 4$+ + - CoA'^-) - - h^-] + 

d^xlc^ix) + JxX + 9'^0[x]\ (2.9) 
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where C describes a contour in the complex time plane as follows: along the forward branch {ti, +00) the fields 
and sources are $+, x"*", along the backward branch (00, U) the fields and sources are <i>~, x^, and along the 
Euclidean branch [ti.ti — if}) the fields and sources are $ = 0;x'',</^- Along the Euclidean branch the interaction 
term vanishes since the initial density matrix for the field x is assumed to be that of thermal equilibrium. The contour 
is depicted in fig. ^ 



ti 



X" 



T+ T+ 



00 



X ' '^x 

FIG. 1: Contour in time for the non-equilibrium path integral representation. 



The linear term /i<i>^ is a counterterm that will be required to cancel the linear terms (tadpole) in $^ in the non- 
equilibrium effective action. This issue will be discussed below when we obtain the non-equilibrium effective action 
for the field $ after integrating out the field(s) x- 

The trace over the degrees of freedom of the x field with the initial equilibrium density matrix, entail periodic 
(for bosons) or antiperiodic (for fermions) boundary conditions for x along the contour C. However, the boundary 
conditions on the path integrals for the field <& are given by 



and 



<i?+(x, i = 00) = <1> (x,t = oo) 



(2.10) 



(2.11) 



The reason for the different path integrations is that whereas the x field is traced over with an initial thermal 
density matrix (since it is taken as the "bath"), the initial density matrix for the $ field will be specified later as part 
of the initial value problem. The path integral over x leads to the infiuence functional for $^[l3- 

A. Tracing over the "bath" degrees of freedom 

As far as the path integrals over the bath degrees of freedom x is concerned the fields $^ are simply c-number 
sources. The contour path integral 



(2.12) 



is the generating functional of correlation functions of the field x in presence of external c-number sources $^ (the 
sources , 
namely 



sources generate the correlation functions via functional derivatives and are set to zero at the end of the calculation), 



/ 



(2.13) 
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Note that the expectation value in the right hand side of eqn. (|2.13() is in the equihbrium density matrix of the 
field X- The path integral can be carried out in perturbation theory and the result exponentiated to yield the effective 
action as follows 



/ .s/,d^.*OM\ I d4^<i,(^)/0[^](^A I ^4^ I dV$(x)$(a:')(0[x]W0M(x')) +0{9^) (2.14) 

X •/ f X ^ 1 1 X 



Ml 

X '2 Jc Jc 'X 

This the usual expansion of the exponential of the connected correlation functions, where this series is identified 
with 



X 



(2.15) 



and where the influence functional\l^ Li/[$+, $ ] is given by the following expression 

L,f[<f+,^-]=g J^d''xMx){0[xKx))x+i^ J d^x j d^x'<^{x)^{x'){0[x]{x)0[x]{x'))^,eon + 0{g'') (2.16) 



In detail, the integrals along the contour C stand for the following: 



^d''x'^{x){0[x]{x))^^ Jd'xJ^ dt[^+ix,t){0[x+Kx))^~<p-ix,t){0[x-]{x))^] (2.17) 

r r poo p poo 

^d'x y^dV<J>(x)<I>(x')(OM(^)0[x](^'))x,con= J d^^ dt J d^^' dt'[^+{x)<S>+{x'){0[x+]{x)0[x+]{x'))x,cou 

+ <p-{x)<^>-{x')(0[x-]{x)0[x-]{x'))^^con 

- ^+{x)^-{x')(0[x+]{x)0[x-]{x'))^^con 

- ^-{x)<^>+{x'){0[x-]i^)0[x+]{x'))^,con\ (2.18) 

Since the expectation values above are computed in a thermal equilibrium translational invariant density matrix, 
it is convenient to introduce the spatial Fourier transform of the composite operator O in a spatial volume V as 

O-.it) = ^7 / d'xe^'-'0[x{x,t)] (2.19) 



in terms of which we obtain following the correlation functions 

{Oj:it)) ^ {0+{t)) (0-{t)) = Tre-P"^ 0-^{t) (2.20) 

{Oj:{t)0_^^{t')) = {0.{t)0+_.{t'))=TTO_j:{t')e-P"^0^^{t)=g>{t-t')=g-+{t,t') (2.21) 

(0_,r(i')0£(i)> - (Ot(t)0:,-(t')> =TrO,~(t)e-^«-^0_,,(i') =a,<(t-i') (2-22) 

{TO^^{t)o_j:{t')) = g>{t - t')Q{t - 1') + g<{t t')Q{t' -t) = g++{t,t') (2.23) 

{fOf:it)o_^,{t')) ^g>it-t')eit' -t) + g<{t-t')e{t-t') = g--{t,t') (2.24) 

The time evolution of the operators is determined by the Heisenberg picture of H^, namely Oj:{t) — 

g«^fx(*~*i)(9-.(i^)g-*^fx(*-*i) Because the density matrix for the bath is in equilibrium, the correlation functions 
above are solely functions of the time difference. These correlation functions are computed exactly to all orders in 
the couplings of the bath fields amongst themselves. 

These correlation functions are not independent, but obey 

Qt^M + Qrii^t') - QI^M - Qt~M = (2.25) 
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The non-equilibrium effective action is given by 



dtd^x [£o.$(*+) + h<^+ - £o,*(*") - h<!>-] + L.,/[$+, (2.26) 

where we have set the sources for the fields <i>^ to zero. 
The choice of counterterm 

h = -{0{x,t)) (2.27) 

cancels the terms linear in $^ (tadpole) in the non-equilibrium effective action. 

In what follows we take ti = without loss of generality since (i) for t > ti the total Hamiltonian is time independent 
and the correlations will be solely functions oit — ti, and (ii) we will be ultimately interested in the limit t ^ ti when 
all transient phenomena has relaxed. In terms of the spatial Fourier transform of the fields defined as in eqn. 
we find 



nOO 

$+(t)C?++(i,t')<f>+^-(t') + *£(i)efe--(i,<')$:,~(^') 



A; 

dt / dt' 



_9_ 
2 Jo 

-*^(i)5.t-(t,t')'i>:£(i') - 'f^Jt)g-+{t,t')<P+-it')] I (2.28) 

where all the time integrations are in the interval < t < oo. 

A similar program has been used recently to study the relaxation of scalar fields[23| as well as the photon production 
from a quark gluon plasma in thermal equilibrium jl^]. 

B. Stochastic description: generalized Langevin equation. 

As it will become clear below, it is more convenient to introduce the Wigner center of mass and relative variables 
*(x,i) = i($+(.x,i) + $-(f,i)) ; R{x,t)={^+ix,t)-<^-ix,t)) (2.29) 
and the Wigner transform of the initial density matrix for the $ field 

(2.30) 

The boundary conditions on the $ path integral given by (|2.11f) translate into the following boundary conditions on 
the center of mass and relative variables 

*(f , t = 0) = *j ; R{x, t = 0) = R, (2.31) 
furthermore, the boundary condition (|2.10|) yields the following boundary condition for the relative field 

R{x,t ^ oo) =0. (2.32) 

This observation will be important in the steps that follow. In terms of the spatial Fourier transforms of the 
center of mass and relative variables H2.29|l introduced above, integrating by parts and accounting for the boundary 
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conditions H2.31|l . the non-equilibrium effective action H2.28|l becomes: 

iL,ff[^,R] = J d^^{-^i^_£(#,,(^) + (fc' + m2)^'fc(^))} 

k 

- dtj^ dt' |ii?_,-(i)i?,^(t')/Cfc(i - t') + R_S)^^k(t - 
+ / d^xRi{x)^{x,t = 0) 



(2.33) 

where the last term arises after the integration by parts in time, using the boundary conditions (|2.31f) and H2.32|l . 
The kernels in the above effective Lagrangian are given by (see eqns. (|2. 2 112. 2411 ') 



/Cfc(i-t') 



[Gl{t~t') + g<{t-t')] 



i^^{t~t') = g^[g>{t-t')^g<{t-t')\Q{t-t') = i^k{t~t')Q{t~t') 
The term quadratic in the relative variable R can be written in terms of a stochastic noise as 

exp[-^jdtj dt'R_^{t)ICk{t~t')R^{t')} = J I?eexp{-i J dt J dt' ^^W^k' ^ t')^-kit') 

+t J dt e_£(i)i?fe(t)} 

The non-equilibrium generating functional can now be written in the following form 

Z = j D^, j DU,J V-^VRV^ >V(«',;n,)i:»i?,e*^'''^-"'(^)(n'(^)-*(^'*=°»-p[e] 



(2.34) 
(2.35) 



exp 
exp 



dt Rt{t) 



dt 



oc 



*,-(t) + (fc2 + m')'ff:{t) + / dt' J:j:it - t')*fc(t') - m 



dt' ^^^{t)K^\t~t')(,_~^{t') 



(2.36) 



(2.37) 



(2.38) 



The functional integral over Ri can now be done, resulting in a functional delta function, that fixes the boundary 
condition '^{x^t = Q) = Iii{x). 

Finally the path integral over the relative variable can be performed, leading to a functional delta function and the 
final form of the generating functional given by 



^M) + (fc2 + m2)*-(t) + / dt' Efe(t - t')*fc(<') - CfeW 



(2.39) 



(2.40) 



with the boundary conditions on the path integral on ^ given by 

5-(f,t = 0) = *,(f) ; i'{x,t^Q)^U,(x) 

where we have used the definition of S^(i — t') in terms of Efe(t — t') given in equation 

The meaning of the above generating functional is the following: in order to obtain correlation functions of the 
center of mass Wigner variable 5" we must first find the solution of the classical stochastic Langevin equation of motion 



(2.41) 



^j:{t) + (fc2 + m')^j:{t) + / dt' Efe(t - <')*,-(*') = 

Jo 

*£(^-0)^v|/^.j; vi;^(t^O)=m£ 



for arbitrary noise term ^ and then average the products of \1/ over the stochastic noise with the Gaussian probability 
distribution V[£] given by (|2.38|l . and finally average over the initial configurations '^i{x);IIi{x) weighted by the 
Wigner function W{'^i,Ili), which plays the role of an initial phase space distribution function. 

Calling the solution of H2.41|l 'i'j:{t;£^;'i>i;Ili), the two point correlation hmction, for example, is given by 

{^_^{t)^j:{t')) = J V[^]V[^] J D^, J DU, W(vI/,;n,)^fc(t;^;^,;nO*_,-(i';^;^.;nO (2.42) 
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We note that in computing the averages and using the functional delta function to constrain the configurations of 
^ to the solutions of the Langevin equation, there is the Jacobian of the operator /dt^ + (fc^ + to^) + / dt'Y7f^^ {t — t') 
which however, is independent of the field and cancels between numerator and denominator in the averages. 

This formulation establishes the connection with a stochastic problem and is similar to the Martin-Siggia-Rose!22| 
path integral formulation for stochastic phenomena. There are two different averages: 



The average over the stochastic noise term, which up to this order is Gaussian. We denote the average of a 
functional JF[^] over the noise with the probability distribution function P[^] given by eqn. H2.38|l as 

(t^ia» . iyM- (2.43) 

Since the noise probability distribution function is Gaussian the only necessary correlation functions for the 
noise are given by 

((CsW)) =0, {{i-^(t)ij:,{t'))) ^K.^(t-t')5\k + k') (2.44) 

and the higher order correlation functions are obtained from Wick's theorem. Because the noise kernel JCk{t — 
t') ^ S{t — t') the noise is colored. 

The average over the initial conditions with the Wigner distribution function 'W{^i,Ili) which we denote as 



In what follows we will consider a Gaussian initial Wigner distribution function with vanishing mean values of 
Hi with the following averages: 

1 

2Wk 



= [1 + 2AA^fc] (2.46) 



Wk 



nan,-£-^[l + 2AA,,,] (2.47) 



na*.-fe + *an.-*;-o (2.48) 

where Wk is a reference frequency. Both Wk and JVb^k characterize the initial gaussian density matrix. Such 
a density matrix describes a free field theory of particles with frequencies Wk- The averages (|2.46I2.47|) are 
precisely the expectation values obtained in a free field Fock state with JVb^k number of free field quanta of 
momentum k and frequency Wk or a free field density matrix which is diagonal in the Fock representation of a 
free field with frequency Wk- This can be seem simply by writing the field and canonical momentum in terms 
of the usual creation and annihilation operators of Fock quanta of momentum k and frequency Wk - While this 
is a particular choice of initial state, we will see below that the distribution function becomes insensitive to it 
after a time scale longer than the quasiparticle relaxation time. 

The average in the time evolved full density matrix is therefore defined by 



(•••)^ ((•••)) • (2.49) 

C. Fluctuation and Dissipation: 

From the expression (|2.35|) for the self-energy and the Wightmann functions H2.21I2.22|I which are obtained as 
averages in the equilibrium density matrix of the x fields (bath), we now obtain a dispersive representation for the 
kernels JCk{t — t'); T,j^{t — t')- This is achieved by explicitly writing the expectation value in terms of energy eigenstates 
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of the bath, introducing the identity in this basis, and using the time evolution of the Heisenberg field operators to 
obtain 

/oo 
d.;a|He-(*-*') (2.50) 
-OO 



g'G^it-t') = I a5(c.) e-(*-*') (2.51) 



with the spectral functions 



a|(a>) = ^^e-'5^"(n|O,,(0)|m)(m|O_£(0)|n)(5(co- (£;„-^„0) (2-52) 



b 



= ^y,e-^''^-{n\0_j:mm){m\Oj:mn) 5{ij - {E„, - En)) (2.53) 

m,n 

where Zt, — Tr e^^^^ is the equilibrium partition function of the "bath" . Upon relabelling m <-> n in the sum in the 
definition we find the KMS relation 23, 24J 

a<{L.)=a>{-uj)^e^-<j>{i.) (2.54) 

where we have used parity and rotational invariance in the second line above to assume that the spectral functions 
only depend of the absolute value of the momentum. 

Using the spectral representation of the Q{t — t') we find the following representation for the retarded self-energy 



J-oo 27r 



(2.55) 



with 



£-(Mo) = r ^..1^^^^^^^^ (2.56) 
J_oo w - fco + 

Using the condition (|2.54|l the above spectral representation can be written in a more useful manner as 



S«(A:,fco) = — / du: , (2.57) 

where the imaginary part of the self-energy is given by 

ImE^(fc,w) = ^cT>(cj) [e^"" - l] (2.58) 

and is clearly positive for ui > 0. Equation (|2.54|) entails that the imaginary part of the retarded self-energy is an odd 
function of frequency, namely 

ImS^(fc, u) = -ImE^(fc, -uj) . (2.59) 
The relation (|2.58() leads to the following results which will be useful later 

ar:{Lu) = ilmS^(fc,w)n(w) (2.60) 

TT 

a<{uj) = ilmS^(fc,w) [1 -hn(w)] (2.61) 

TT 
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k,k, 



,I\0 




£(k,ko) 



FIG. 2: Self-energy of $ to lowest order in but to all orders in the couplings of the fields x amongst themselves. The external 
fines correspond to the field $. 

Similarly from the definitions H2.34|l and (|2.50I2.51|I and the condition H2.54|) we find 

^k{t-t') = r ^e*'=o(*-*')^(fc,fc^) (2.62) 
J-oo 27r 

^(fc,fco) = ^c7>(fco) [e'^'=« + l] (2.63) 
whereupon using the condition H2.54|l leads to the foUowint generalized form of the fluctuation-dissipation relation 



/C(fc, ko) = ImE'"(fc, fco) coth 



(2.64) 



Thus we see that ImE^(fc, ka) ; /C(fc, fco) are odd and even functions of frequency respectively. 
For further analysis below we will also need the following representation (see eqn. 



/OO -POO 
^^u(t-t') ^^>(^^ _ ^<(^)j duj^- e*"(*-*')lml]^(fc,w)dw (2.65) 
-OO ^ J — OO 

whose Laplace transform is given by 



S(fc,s) = r dte-^'l^kit) = -- r ^^^^^duo (2.66) 

Jo J-oo ^ + «S 

This spectral representation, combined with 12.57|l lead to the relation 

E^(fc,fco) = S(fc,s ifco + e) (2.67) 

We highlight that the self-energy Y,^{k^ fco) as well as the fluctuation kernel /C(fc, fcg) are to all orders in the couplings 
amongst the fields x but to lowest order, namely 0{g^) in the coupling between the field $ and the fields The 
self-energy is depicted in fig.Q. 

D. The solution: 

The solution of the Langevin equation (|2.41|l can be found by Laplace transform. Defining the Laplace transforms 

^ /-OO 

^j:{s) ^ \ dte-^'^^,{t) (2.68) 

^ /"OO 

4^(s)EE / dte-''i^{£) (2.69) 
Jo 

along with the Laplace transform of the self-energy given by eqn. (|2.66(l we find the solution 

,T, ^ \ i,k ' i.k ' ■> 2 7 2,2 /o ■7n^ 

^fcls) = — ^ — — I ^k^k + TO (2.70) 
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where we have used the initial conditions H2.40|l . The solution in real time can be written in a more compact manner 
as follows. Introduce the function fk{t) that obeys the following equation of motion and initial conditions 



ft 

,2 



fk{t)+u;lfk{t)+ I dt' ^k{t-t')fk{t')^Q ; /(t = 0) = 0; /^(t - 0) = 1 (2.71) 
whose Laplace transform is given by 



^^^'^ = ^ (2.72) 

In terms of this auxiliary function the solution of the Langevin equation (|2.41|) in real time is given by 



vl/fc(i;vl/,;n,;e)=vl'^£ A.W+n^,-/fc(i)+ / fk{t ~ t') ^^{t')dt' (2.73) 

JO 

For the study of the number operator below we will also need the time derivative of the solution, given by 

i>k{t; n,; ^ = £ fk{t) + fk{t) + f fk{t - t') ^^{t')dt' (2.74) 

where we have used the initial conditions given in eqn. (|2.73(l . From eqn. I|2.7()(l it is clear that the solution (|2.73(l 
represents a Dyson resummation of the perturbative expansion. 

The real time solution for f{t) is found by the inverse Laplace transform 

/fc(t) = / (2.75) 

where C stands for the Bromwich contour, parallel to the imaginary axis in the complex s plane to the right of all the 
singularities of /(s) and along the semicircle at infinity for Res < 0. The singularities of /(s) in the physical sheet 
are isolated single particle poles and multiparticle cuts along the imaginary axis. Thus the contour can be deformed 
to run parallel to the imaginary axis with a small positive real part with s = ia; + e;— oo<u;<oo, returning parallel 
to the imaginary axis with s ^ iuj — e; oo > uj > —oo, with e = 0+ as depicted in fig. 




FIG. 3: General structure of the self-energy in the complex s-plane. The dashed regions correspond to multiparticle cuts namely 
ImE^(fc, s = io; + e) 7^ 0. The dots depict isolated poles. 
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From the spectral representations p. 5812. 66(1 1 one finds that S(A;, s = iuizte) = ReE^(fc, uj) ± ImE^(fc, uj) and using 
that lmS^{k,uj) = —lmY,^{k, —uj) we find the following solution in real time 



fk{t)= I s\n{ujt) p{k,uj;T)dLjj , 
where we have introduced the spectral density 



(2.76) 



p{k,uj;T) = - 

TT 



lm^^{k,uj;T) + 2ue 



^ 1 2 r n 2 ' 

-ujl-Yiei:"{k,uj]T) + ImE«(fc,tj;r) + 2we 



(2.77) 



and we have made explicit the temperature dependence of the self-energy. 

We have kept the infinitesimal 2ti;e with e ^ 0^ since if there are isolated single particle poles away from the 
multiparticle cuts for which ImE^(A;,s) = then this term ensures that the isolated pole contribution is accounted 
for, namely 



1 

TT 



2uje 



Lu^-Lol- ReY.^{k, uj) + [2ujef 
The initial condition fk{t — 0) — ^ leads to the following sum rule 



sign(w) S uj^ -ujI- ReE^(fc, uj) 



(2.78) 



duj 

TT 



Iml]^(fc,w) + 2uje 



ReE«(fc,c 



= 1 



ImE^(fc,w) + 2uie 



(2.79) 



III. COUNTING PARTICLES: THE NUMBER OPERATOR 



In an interacting theory the definition of a particle number requires careful consideration. To begin with, a 
distinction must be made between physical particles that appear in asymptotic states and can be counted by a 
detector, from unstable particles or resonances which have a finite lifetime and decay into other particles. Resonances 
are not asymptotic states, do not correspond to eigenstates of a Hamiltonian and their presence is inferred from virtual 
contributions to cross sections. In an interacting theory virtual processes turn a bare particle into a physical particle 
by dressing the bare particle with a cloud of virtual excitations. Physical particles correspond to asymptotic states 
and are eigenstates of the full (interacting) Hamiltonian with the physical mass. These physical particles correspond 
to real poles in the Green's functions or propagators in the complex frequency plane. In the exact vacuum state, 
the propagator of the field associated with the physical particles features poles below the multiparticle continuum at 
the exact frequencies and with a residue given by the wave function renormalization constant Z. The wave function 
renormalization determines the overlap between the bare and interacting single particle states. Lorcntz invariance of 
the vacuum state entails that the exact frequencies are of the form ilk — + nip, where mp is the physical mass 
and that the wave function renormalization is independent of the momentum k. In asymptotic theory, the spatial 
Fourier transform of the field operator ^^it) obeys the (weak) asymptotic condition 



/J^""«L|0> . yje--!!,), (3.1) 

where jl^:) is the state with one physical particle. 

In a medium at finite temperature there are no asymptotic states, each particle, even when stable in vacuum 
acquires a width in a medium either by coUisional processes (coUisional broadening) or other processes such as Landau 
damping. The width acquired by a physical particle in a medium is a consequence of the interaction between the 
physical particle and the excitations in the medium. In particular the medium-induced width is necessary to ensure 
that physical particles relax to a state of thermal equilibrium with the medium. The relaxation rate is a measure of 
the width of the particle in the medium. Therefore in a medium a physical particle becomes a quasiparticle with a 
medium modification of the dispersion relation and a width. 
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Thus the question arises as to what particles are "counted" by a definition of a distribution function, namely, a 
decision must be made to count either physical particles or quasiparticles. 

One can envisage counting physical particles by introducing a detector in the medium. Such detector must be 
calibrated so as to "click" every time it finds a particle with given characteristics. A detector that has been calibrated 
to measure physical particles in a scattering experiment for example, will measure the energy and the momentum 
(and any other good quantum numbers) of a particle. Every time that the detector measures a momentum k and an 
energy D,k determined by the dispersion relation of the physical particle (as well as other available quantum numbers), 
it counts this "hit" as one particle. 

Once this detector has been calibrated in this manner, for example by carrying out a scattering experiment in the 
vacuum, we can insert this detector in a medium and let it count the physical particles in the medium. 

Counting quasiparticles entails a different calibration of the detector which must account for the properties of the 
medium in the definition of a quasiparticle. The first obstacle in such calibration is the fact that a quasiparticle does 
not have a definite dispersion relation because its spectral density features a width, namely a quasiparticle is not 
associated with a sharp energy but with a continuum distribution of energies. How much of this distribution will be 
accepted by the detector in its definition of a quasiparticle, will depend on the filtering process involved in accepting 
a quasiparticle, and so cannot be unique. Therefore statements about measuring a distribution of quasiparticles are 
somewhat ambiguous. 

In this article we focus on the first strategy, by counting only physical particles. Hence we propose a number 
operator that "counts" the physical particle states of mass mp that a detector will measure for example in a scattering 
experiment at asymptotically long times. Asymptotic theory and the usual reduction formula suggest the following 
definition of an interpolating number operator that counts the number of physical (stable) particles in a state 

Mt) = {hj:{t)h_^{t) + r!2$-(i)$_~(t)} - Ck (3.2) 

where Z is the wave function renormalization, namely the residue of the single (physical) particle pole in the exact 
propagator, Qk = \/fc2 + rrip is the renormalized physical frequency and the normal ordering constant Ck will be 
adjusted so as to include renormalization effects. In free field theory il.k = ujk — VW + rn^ , Z — 1 , Ck = 1/2. 
However, in asymptotic theory the field $ creates a single particle state of momentum k and mass mp with amplitude 
^/Z out of the exact vacuum. 

The quantity Ck arises from the necessity of redefining the normal ordering for the correct identification of the 
particle number in an interacting field theory. It will be fixed below by requiring that the expectation value of Nk {t) 
vanishes in the exact vacuum state at asymptotically long time. Alternatively this constant can be extracted from 
the equal time limit of the operator product expansion. 

The approach that we follow is to consider an initial factorized density matrix corresponding to a tensor product 
of a density matrix of the field and a thermal bath of the fields x- This initial state will evolve in time with the 
full interacting Hamiltonian, leading to transient phenomena which results in the dressing of the bare particles by the 
virtual excitations. At asymptotically long times the bare particle is fully dressed into the physical particle, and at 
finite temperature, a quasiparticle. The time evolution of the interpolating number operator will reflect this transient 
stage and the dynamics of the dressing of the bare into the physical state. Since the thermal bath is stationary, the 
distribution of physical particles in the bath will be extracted from the asymptotic long time limit of the expectation 
value of the interpolating Heisenbcrg number operator Nk{t) in the initial state. 

The expectation value of Nk (t) is related to the real-time correlation functions of the field <& as follows 

- Ck (3.3) 

t=t' 

where the non-equilibrium correlation functions are given by 

t')+g<it,t')eit'-t) (3.4) 
t)+g<{t,t')Qit-t') (3.5) 

(3.6) 
(3.7) 



1 

m-z 



d_d_ 

dtdt' 



9l{t,t')+g<{t,t') 



{n(^W-k^t'))^9t{txm- 

($,(i)$:-(0)=.9fc>(t,0e(t'- 

{>^lJt'm{t))^g<{t,t') 



The symmetrized definition of the expectation value (|3.3|l has been chosen for convenience because, as it is shown 
in detail below, it is related to a simple correlation function of the center of mass field 4'fc(i) given by eqn. H2.29|l . 
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However we could have chosen any other definition of the equal time correlation function, such as (t,t); g^(t,t) 
or any combination thereof that does not involve a Heaviside step function in time for which the time derivatives 
will yield spurious delta functions. It is a straightforward exercise with the density matrix to show that all of these 
alternative definitions are equivalent in the equal time limit, since these do not involve discontinuous functions of 
time. 

In terms of the center of mass field ^fe(i) — i^^{t) + ^^{t))/'^ introduced above it is straightforward to find that 
the correlation function in the bracket in (|3.3() is given by 



(%(i)*_£(i')> = ^k>(t,0+3fc<(i,t')] 

and the occupation number can be written in terms of the center of mass Wigner variable introduced in eqn. 
as follows 



(3.8) 



(Nkit)) 



1 



(3.9) 



where the expectation values are obtained as in eqn. (|2.49() and ^{t) is the solution of the Langcvin equation given 
by J2.V3l2.V4L 

A straightforward calculation implementing eqn. H2.49|l writing the noise in terms of its temporal Fourier transform 
and using the Fourier representation of the noise kernel (|2.62|) leads to the following result 



Nk(t) = {Nk{t)) = 



1 



where we have introduced 



2nkZ ]2Wk 



1 



[1 + 2Ub,k] 



+ 



m + {ni + wi) m + niwi m 



(3.10) 



nk{uO,t) = / dTfk{T)e 



Jo 



(3.11) 
(3.12) 



Jk{t) is given in eqn. (|2.76() and the fluctuation kernel /C(fc, w) is given by eqn. H2.64|l . 

The result H3.1()(l for the time evolution of the distribution function, along with the expressions 1)3. 1113. 12|l clearly 
highlights the non-Markovian nature of the evolution. The integrals in time in (|3.11I3.12|I include memory of the 
past evolution. This is one of the most important aspects that distinguishes the quantum kinetic approach from the 
usual Boltzmann equation. We will contrast these aspects in section (jVll . 



A. Counting physical particles in a thermal bath. 

In the vacuum the spectral density of the field $ which describes a physical particle is depicted in fig. (0J. It 
features isolated poles along the real axis in the physical sheet in the complex frequency (w) plane at the position of 
the exact single particle dispersion relation f2fc with < \uith\ where uth is the lowest multiparticle threshold. 

As mentioned above, in a medium stable physical particles acquire a width as a consequence of the interactions 
with physical excitations, and become quasiparticles. The width can originate in several different processes such as 
collisions or Landau damping. The poles move off the physical sheet into the second (or higher) Riemann sheet in 
the complex lu plane, thus becoming a resonance. This is the statement that there are no asymptotic states in the 
medium. 

The analytic structure of the spectral density at finite temperature is in general fairly complicated. While at zero 
temperature the multiparticle thresholds are above the light cone |a;| > k, at finite temperature (or density) there 
appear branch cuts with support below the light cone p3. 124 . 125ll26| . However a general statement in a medium is that 
the poles associated with stable particles in vacuum (along the real axis in the physical sheet) move off the physical 
sheet and the spectral density does not feature isolated poles but only branch cut singularities in the physical sheet, 
associated with multiparticle processes in the medium. 
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FIG. 4: Spectral density pk{ijJ,T = 0) for stable particles. The dots represent the isolated poles at ±57^ and the shaded regions 
the multiparticle cuts, ujth is the lowest multiparticle threshold. 

In perturbation theory the resonance is very close to the real axis (but in the second or higher Riemann sheet) and 
the width is very small as compared with the position of the resonance. We will study a particular example in the 
next section. 

In perturbation theory the spectral density p(k, uj, T) (12.77(1 features a sharp peak at the position of the quasiparticle 
"pole" which is determined by 

Wl{T) ReS«(fc, Wfc(r); T) = (3.13) 

Near the quasiparticle "poles" the spectral density is well described by the Breit-Wigner approximation 



PBw{k,uj;T) 



sign(c^)rfe(r) 



Zk{T) 1 

2Wk{T) 7r(|c^|->Vfc(r))2 + r2(r) 



where Wk{T) is determined by eqn. H3.13|) and the finite temperature residue and width are given by 



(3.14) 



1 - 



1 



aReI]^(fc,w;T) 



2Wk{T) 



duj 



(3.15) 



rfc(T) 



ImS]^(fc,>Vfc(r);r) 
2Wfe(T) 



(3.16) 



At zero temperature of the bath, the (quasi) particle dispersion relation Wk{T) is identified with the dispersion 
relation of the stable physical particle, namely the "on-shell" pole, the residue 2k{T) is identified with the wavefunction 
renormalization constant Z which is the residue at the on-shell pole for the physical particle, and the width vanishes 
at zero temperature since the particle is stable in the vacuum, namely 



Wfe(r = 


0) 


= n 




0) 


^ z 




0) 


= 



(3.17) 
(3.18) 
(3.19) 



In the Breit-Wigner approximation the real time solution is easily found to be 



frit) ^ Z,{Tf-^^^^^e^^^(^^^ (3.20) 

This solution describes the relaxation of single quasiparticles, where Wfe(i) is the quasiparticle dispersion relation 
and Tk{T) is the quasiparticle decay rate. 
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The asymptotic long time limit of the distribution function (|3.10|l is obtained by using the following identities 



oo 



7^fe(c^,oo)= / e-^("— = f,{s = tu; + e) (3.21) 

JQ 

poo 

^fc(t^,oo) = / e-'^^"''^'jk{t)dt = iiofk{s ^loj + t) (3.22) 



where fk{s) is the Laplace transform of fk{t) given by eqn. H2.72|l and in 13.22|l we have integrated by parts, used 
the initial condition /fc(0) — and introduced a convergence factor e ^ 0+. Hence the expectation value of the 
interpolating number operator in the asymptotic long-time limit is given by 

/■oo / 2 I o2 \ 

iVfe(oo)=y^ f 1 [l + 2n(u;)]p(fe,a;,T)(ia; -Cfc , (3.23) 

where n{u)) is the Bose-Einstein distribution function and we have used the fluctuation-dissipation relation 12.64|l as 
well as eqn. (|2.72l) which lead to the p{k,uj,T) in (|3.23() . The dependence of the asymptotic distribution function on 
the spectral density is a consequence of the fluctuation-dissipation relation (|2.64(l as well as the non-Markovian time 
evolution as displayed in H3.21I3.22|) . 

The real time solution H3.20|l clearly reveals that the asymptotic limit is reached for t > Tk ~ l/2rfc(r) where Tk{T) 
is the quasiparticle relaxation rate. The distribution function at t >> does not depend on the initial distribution 
Mb,k or the reference frequencies Wk- Therefore at times longer than the quasiparticle relaxation time the distribution 
function becomes independent of the initial conditions. This is to be expected if the state reaches thermal equilibrium 
with the bath, since in thermal equilibrium there is no memory of the initial conditions or correlations. 

The integral term in the asymptotic distribution H3.23|l is easily understood as full thermalization from the following 
argument. 

Let us consider the correlations functions g^{t,t');g^{t,t') given by eqns. H3.6I3.7|) . In thermal equilibrium they 
have the spectral representation 

g>(i,0 = I p>{k,u;;T)e^-'^'-'">duj (3.24) 



9k 



<{t,t') = / p<(fc,w;T)e^"(*-*')dw (3.25) 



where 



p>ik,uj;T) = -^y^e-^^-{n\^^{0)\m){m\^_^iO)\n)S{io^iE.r,^Em)) (3.26) 
p<(fc,c^;T) - ^Ve-''^'"(n|$_j(0)|m)(m|%(0)|n)5(cc> -(£;„- £;„)). (3.27) 

ra.n 

Where Zt is the thermal equihbrium partition function. A straightforward re-labelhng of indices leads to the 
relation 



p<(k,u:-T)^p>(k-uj-T)^ e^^p~> {k, to; T) (3.28) 

The spectral density is given by 

p(fc, Lo- T) = p< (fc, T) - p> (fc, Lo- T) (3.29) 

leading to the relations 



p^{k,uj;T) — p{k,Lu;T) n{uj) 
p<(fc,w;T) = pik,u;;T) [1 + n(a;)] , 



(3.30) 
(3.31) 
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where n{uj) = l/[e^" - 1]. 

Therefore in thermal equilibrium the expectation value of the operator term in eqns. (|3.2I I3.3|l is given by 



t=t' 



k 



[l + 2n{Lu)] p{k,uj;T)duj, 



(3.32) 



which is precisely the integral term in the asymptotic limit given by eqn. (|3.23() . Therefore the expression (|3.23() 
indicates that the excitations of the field $ have reached a state of thermal equilibrium with the bath. The normal 
ordering constant Ck in H3.23|l is a subtraction necessary to redefine normal ordering in the interacting theory and is 
defined from the operator product expansion to yield vanishing number of particles in the vacuum. 

While the asymptotic long time limit can be obtained directly from the spectral representation of the interpolating 
number operator in the equilibrium state, the real time formulation in terms of the non-equilibrium effective action 
has two advantages: i) it makes explicit the connection with the fluctuation dissipation relation and clearly states 
that the equilibrium abundance is determined by the noise correlation function of the bath, ii) the real time dynamics 
clearly shows thermalization on time scales t > Tk- These statements would not be immediately recognized from the 
equilibrium spectral representation. 

The result (|3.23l) becomes more illuminating in the narrow width approximation where the Breit-Wigner approxi- 
mation for the spectral density (|3.14() is supplemented with the narrow width limit [T) which leads to 



p{k,uj;T) 



which in turn leads to the approximate result 



sign{uj)6{\uj\~Wk{T)) 



(3.33) 



Nk{oo) 



Z \2Wk{T)nk 



-n(Wfc(r)) 



(3.34) 



Obviously the zero temperature pole flk and residue Z and their finite temperature counterparts WkiT), Zk{T) 
differ by terms that are of order g^, namely perturbatively small, therefore in the narrow width approximation, which 
itself is a result of the weak coupling assumption one could write 



Nkioo) - n(f7fe) 



(3.35) 



Thus choosing the normal ordering factor Ck — 1/2-1- 0{g^) would lead to the conclusion that the physical particles 
are distributed in the thermal bath with a Bose-Einstein distribution function with the argument being the physical 
pole frequency (at zero temperature). Furthermore the normal ordering constant Ck ~ 1/2 is identified with the usual 
normal ordering of the number operator in the free field vacuum. 

In order to understand in detail the perturbative correction we have to first decide on what are ilk , Z,Ck- The 
importance of the perturbative corrections cannot be underestimated, if the temperature of the bath is much smaller 
than Qk the distribution function n{Qk) ^ 1 and the perturbative corrections can be of the same order or larger. 
What should be clear from the above discussion is that in order to make precise the perturbative correction to the 
abundance, we must specify unambiguously what is being counted. 



1. Physical particles in the vacuum 

The next step is to define flk,Z,Ck- As it was emphasized above, the number operator that we seek counts physical 
particles. These are stable excitations off the full vacuum state of the theory and are associated with isolated single 
particle poles in the spectral density at zero temperature. 

The zero temperature limit of the asymptotic distribution function 13.23|l is 

/■oo I q2 \ 

7Vfc(oo;T = 0) l^-^\p{k,u;,T = 0)du;~Ck, (3.36) 
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At T = the spectral density features the isolated single particle poles away from the multiparticle continuum as 
depicted in fig. 10} . The contribution from the single particle poles to the zero temperature spectral density is given 
by eqn. H2.78|l . therefore we write 

pik, w, T = 0) = sign(tj) d{Lu - nk) + pc{k, u;,T^O), (3.37) 

where pc{k,uj,T = 0) is the continuum contribution with support for \lu\ > ujth, where LUth is the lowest multiparticle 
threshold, and the position of the isolated pole satisfies 

nl-ul-ReT.'^ik^nk) = (3.38) 

At zero temperature Lorentz covariance implies that r2| = fc^ + mp, where mp is the pole mass of the physical 
excitations (asymptotic states). 

The residue at the single (physical) particle pole, Z, is given by 



1 

Z 



^ 1 aReI]«(fc,w;T) 
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(3.39) 



Introducing the zero temperature form of the spectral density (|3.37|1 in the sum rule (j2.79|) the following alternative 
expression is obtained. 

Z = l-2 Lu pc{k,uj,T ^Q)doj (3.40) 
Therefore the asymptotic distribution of particles in the vacuum is given by 

1 , r fio^ + nl^ 



A^fc(oo;T-0) = -+y^ |^___fejp^(fc,,^,r = 0)dL^-Cfe, (3.41) 

The normal ordering term Ck is now fixed by requiring that for T = the vacuum state has vanishing number of 
physical excitations. In other words, by requiring Nk{oo, T = 0) = we are led to 

1 (up' + il'^ \ 

Cfe=o+/ o J'o Pc(fc,c^,r^O)rf^. (3.42) 



2 Jg \ 2Z 

We have kept the lower limit in the integral to be w = for further convenience, however Pc(k^ lu,T ~ 0) vanishes for 

Equations (|3.38|) , H3.39I3.40|I and (|3.42|) determine all of the parameters ilk , Z, Ck for the proper definition of the 
distribution function for physical particles. 

Hence the distribution function of physical excitations in equilibrium with the bath at finite temperature is finally 
given by the simple expression 



AA(fc, T) = Nki^) = I I ""It? M [1 + 2^(^)1 Pik, ^, T) - p4k, c., T - 0) Lc. - i , (3.43) 



2zn 



k 



This is the final form of the asymptotic distribution function of physical particles in equilibrium in the thermal 
bath with flk = \/k^ + ml,; Z;Ck given by equations igiSHl , (E!33) (or (jOUIl . lin^ l respectively. 
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B. Renormalization: 

In renormalizable theories the wavefunction renormahzation constant Z is uhraviolet divergent and the expression 
for the asymptotic distribution function (|3.43|) seems to be ambiguous. However proper renormahzation as described 
below shows that the asymptotic abundance is finite. 

In general the imaginary part of the self-energy can be written as a sum of a zero temperature and a finite 
temperature contribution, the latter vanishing at zero temperature, thus we write 

ImS^(w, fc; T) = ImS^(w, fc) + ImE^(w, fc) (3.44) 

Therefore the real part of the self-energy, which is obtained from the imaginary part by a dispersion relation 
(Kramers-Kronig) can also be written as a sum of a zero temperature plus a finite temperature contribution, 

Re^^iu, k;T)^ V 2k„ , °' ' ' dko = ReS]g^(cj, k) + ReS^(cj, k) (3.45) 

TT Jo k^~ UJ^ 

where V stands for the principal part of the integral, and we have used the fact that lmE^(fco, k\ T) is an odd function 
of fco. Both ImS^(w, k) and ReS^(a;, k) vanish at T = 0. 

The position of the physical pole is obtained at zero temperature from the relation H3.38|l . 

nl -Lol- ReE^(fc, nk) = (3.46) 
The subtracted real part of the self energy is 

ReE^(fc,t^) -ReS^(fc,l]fc) = \l - Z-^[k,uj]\{'^^ - ^l) (3.47) 

where 



■«-a.>)«-Si?) 

This expression for Z^^[k, lo] follows from the zero temperature limit of the dispersive representation in eqn. H3.45|l . 

The function Z[k,Lu = 17^], namely evaluated on the single particle mass shell, is identified with the wave function 
renormalization, or residue at the single particle pole at zero temperature. 

As mentioned above, in renormalizable theories Z[k,uj] is ultraviolet logarithmically divergent, therefore it is con- 
venient to perform yet another subtraction of the integral term in H3.48|l as follows, 

Z-'^[k,Lu]= Z-'^ -na{k,uj), (3.49) 
where Z is the wavefunction renormalization constant, namely the residue at the pole, 



- = l + ip/ 2ko Zt''n^:J dko, (3.50) 



1 p ImI]o^(fco,fc ) 

n^Jo {kl-nlf 

and no(fc,w) is the real part of the twice subtracted self-energy given by 

1 f°° ImS^ffcn k) 

noik,.) = -y-nl)vl 2fco (,g_^,°(^,g°'_;,.). dk, (3.51) 

The two subtractions had been performed on the single particle mass-shell. In a renormalizable theory the integral 
in the twice subtracted real part of the self energy no(A;, w) is ultraviolet finite while the integral in Z~^ is logarith- 
mically divergent. Furthermore the finite temperature parts do not have primitive divergences since all the primitive 
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divergences are those of the zero temperature theory. We emphasize that these expressions are stih functions of the 
bare couphng and any potential divergences arising from coupHng renormaUzation have not yet been accounted for. 
The divergences associated with couphng constant renormahzation wiU be addressed below. 

Combining equations H3.46|l . (|3.47|) and H3.49|l . the spectral density H2.77|l can be written in the following form 



p{k,uj;T) = - 



ImS^(fc,w;r) + 2we 



Z-i(cj2-rj2)_n(fc,cj;r) + ImS«(/fc,w;T) + 2cje 



2 ' 



where 



(3.52) 



n(fc, Lu; T) = _ f7^)no(fc, uj) + Rel]^(w, fc) 
Introducing the renormalized real and imaginary part of the self-energy as 



(3.53) 



nr{k,uj;T) = Zn(fc,w;T) 
ImE^(fc,w;T) = ZImS^(fc,w;T) 

the spectral density l|3.52|l can be written as 



(3.54) 
(3.55) 



where 



p{k,uj;T) = Z pr(k,uj;T) , 



(lu^ -nl)-nr{k,uj;T) + ImE^(fc,a;;r) + 2we 



(3.56) 



(3.57) 



We note that at zero temperature the spectral density Prik, ll!;T — 0) has unit residue at the single physical particle 
pole. 

bmcc both n(fc, Lo; T) an d Iml ] ^(fc,cj ; T) are proportional to , the renormahzation of the real and imaginary part 
of the self-energy in eqns. H3.54|l . (13.551) is tantamount to the renormahzation of the coupling constant^ 



Zg 



(3.58) 



In terms of gr, both Tir{k^ lo\ T) and ImS^(fc, oj; T) are finite since the only counterterms necessary are those of the 
zero temperature theory. Therefore the equilibrium distribution function can be written solely in terms of renormalized 
quantities as follows 



AA(fc,T) = 



2 ill 



[1 -I- 2n(w)] pr{k, UJ, T) - pr,c{k, w, T = 0) > - - 



(3.59) 



This definition of the asymptotic distribution function is one of the main results of this article. 



^ The coupling g in the Lagrangian aheady has the proper renormahzation of the (composite) operator ©[x]- 
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IV. THE MODEL 

The results obtained in the previous section are general and as mentioned above the quantum kinetic effects that 
modify the standard Boltzmann suppression of particle abundance in the medium depend on the particular theory 
under consideration. To highlight the main concepts in a specific scenario, we now consider a theory of three interacting 
real scalar fields with the following Lagrangian density. 



2 2 



2 

E 

i=l 



9^X1 X2 + Ant[XlX2] 



(4.1) 



We will assume that the mutual interaction between the fields xi ? X2 ensures that the fields xi-2 are in thermal 
equilibrium at a temperature T = 1//?. A similar model has been previously studied in ref.|25j| for an analysis of the 
different processes in the medium. 

The particles associated with the field $ will be stable at T = provided mp < Mi + M2, where mp is the 
zero temperature pole mass of the $ particles. In order to study the emergence of a width for the particles of the 
field $ to lowest order in perturbation theory we will consider the case in which Mi > mp + M2 (or alternatively 
M2 > mp + Ml) in this case the quanta of the field xi can decay into those of the field $ and X2- Since the particles 
1, 2 are in a thermal bath in equilibrium the presence of the heavier species (here taken to be that of the field xi) in 
the medium results in a width for the excitations of field <i> through the process of decay of the heavier particle into 
the lighter scalars and its recombination, namely Xi ^ + X2- As will be seen in detail below the kinematics for this 
process is similar to that for Landau damping in the case of massive particles |2^ . 

The relevant quantity is the self-energy of the field <I> which we now obtain to one loop order 0(5^) in the Matsubara 
representation. The one-loop self-energy is given by 



E(z.„,fc) = -.g2 / ^i^G(")(a;„,p)G(«)(c^„-f^„,p + fc) 



(4.2) 



where ujrn, Vn are Bosonic Matsubara frequencies. It is convenient to write the Matsubara propagators in the following 
dispersive form 



4°?(^m,p) 
G(,°)(a;™ + i/„,p + fc) 
Pi(po,p) 

P2{qo,P + k) 



4" 



dpo 
dqo 



1 



Pa - icum ' 

P2{qo,P + k) 
go - iujjn - ivn 
,(1) 



2w^J 



[%o-41-)-%o+^l!^,)], 



2w 



~ p+k 



(2) 



p+k 



p+k' 



pP + Ml; u:'^lj:^\J{p + kf + Ml 



(4.3) 

(4.4) 
(4.5) 

(4.6) 
(4.7) 



This representation allows to carry out the sum over Matsubara frequencies LOm in a rather straightforward 
manner [23, lljl which automatically leads to the following dispersive representation of the self-energy 



T.{k,v.n) = - 

with the imaginary part of the self-energy given by 



duj 



ImS«(fc,tj) 



(4.8) 



ImE-(fc,.).V/70 



dpo I dqo[n{po)-n{qo)]pi{po,p)p2iqo,p + k)S{uj-qo+po) 



(4.9) 
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XhP 



(a) 



X2,P + k X2,P + k 

FIG. 5: Processes contributing to ao{k,Lj),aa{k,ui) (a) and to ai,{k,ui) (b). The inverse processes are not shown. 



where n{q) are the Bose-Einstein distribution functions. From the representation H2.57|l the retarded self-energy 
foUows by analytic continuation, namely 

ko) = S](fc, iy„ =ko- ie) (4.10) 
The imaginary part of the self energy can be written as a sum of several different contributions, namely 



Imi;^(fc, Lu; T) = ao{k, lo) + cra{k, oj; T) + (Ji,{k, T) 
where ao{k^uj) is the zero temperature contribution given by 



(4.11) 



ao{k,uj) 



I (it %) 



p+kj 



P p+k L 

and (Ta(fc, w), tTb(fc, cj) are the finite temperature contributions given by 



(4.12) 



(7a(fc, w;T) 



327r2 ./ ^(.1)^(2), 

P p+k 



(4.13) 



ab{k,ui]T) 



P p+k 



^ p+k' ^ P ' 



\ p p+k) \ p p+k) 



(4.14) 



The processes that contribute to U{j{k, uj) and <Ta{k^ lo) are ^ X1X2 while the processes that contribute to ab{k, uj) 
are xi,2 ^ ^ X2..1 depicted schematically in fig. © 

The details of the calculation of the different contributions are relegated to the appendix. The result is summarized 
as follows: 



cro(fc, w) 



167rQ2 



sign(c^)e[Q^-(Mi+M2)^] {Qy - 2Q\Ml + M^) + {Mt - M^Y ; Q' ^ uj' - k' (4.15) 



We have explicitly displayed the fact that the zero temperature contribution to the imaginary part is manifestly 
Lorentz invariant and solely a function of the invariant mass Q' = uP' ~ k'^ . The finite temperature contributions are 



aaik,Lu;T) 



9 



\QTTk(3 



sign(a;)e[Q2_(Mi + M2)2] 



In 



'l-e-^<' 



Ml ^ M2 



(4.16) 



ab(k,u;T) 



9 



sign(w)e[(Mi-M2)2~Q2i 



IQnkP 



/l_e-'3Kl\ 
\ 1 - e-'^Kl / 



(4.17) 



25 



where 



(4.18) 



The real part of the self energy is obtained from the dispersive form (|2.57|l and can be separated into a zero 
temperature and a finite temperature part as follows 



ReS^(fc, w; T) = ReS^(fc, u) + ReS^(fc, T) 



(4.19) 



with 



ReS]^(fc,w) = 
ReS^(fc,cj;r) = 



_ 1^ r aa{ko,k;T) + at{ko,k;T) 



(4.20) 
(4.21) 



where 7^ stands for the principal part. We note that both ao{k,u!) and aa{k,uj) feature the standard two particle 
threshold above the light cone at the invariant mass = (Mi + M2)^ whereas the finite temperature contribution 
ab{k,u}) has support for invariant mass < (Mi — M2)^ even below the light cone and vanishes at T = 0. In 
the case of massless particles in the loop this contribution is below the light cone and is identified with Landau 
damping j23l 123 . . In particular at zero temperature the isolated poles are at = mp, hence if m% < {Mi — M2)^ 
the physical particle pole is embedded in the multiparticle continuum and moves off the real axis onto the second (or 
higher) Riemann sheet in the complex frequency plane. Because of this the physical particle acquires a width. The 
spectral density for the case rrip < (Mi — M2)^ is depicted in fig. © 



(CO 


















mm. v/m 









FIG. 6: Spectral density p{k,uj,T) for mp < (Mi - M2 



into an unphysical sheet. 



and ujth2 = ^/W+{^hTM2 



^. The shaded areas are the multiparticle cuts with thresholds 
The single particle poles at n| = fc^ + m% moved off the real axis 



1. Zero temperature: 0,k;Z;Ck: 

Using that (To(fco,^) is odd in ko and that it is solely a function of the invariant — kg ~ k"^ for fco > 0, it is 
straightforward to find the following manifestly Lorentz invariant result 



1 



ReT.^{k,u) = V 

./(Mi+Af2)2 



P^-Q^ 



dp^ 



g2 = _ fc2 



(4.22) 
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where we have exphcitly exhibited the two particle threshold in the lower limit. Lorentz invariance requires that the 
single particle pole features the dispersion relation fl^ — k"^ + m?p, and so the equation that determines the single 
particle physical poles, namely eqn. H3.38|l is given by 

m%-m^ - ReS^(Q2 ^ ^2^) = q (4.23) 
From the results of the previous section (see eqn. H3.50|l ) the wave function renormalization constant is given by 



z- 



(4.24) 



Separating the residue at the physical particle pole and following the steps described in section Ijlll Bll the renor- 
malized spectral density H3.56I3.57|I at zero temperature can now be written in the following simple form 



Prik,uj;T = 0) 



1 



iQ^-mj,)(l-UoAQ 



[(7o,.(Q2) + 2u;e]^ 



= sign(u;) S{uj'^-nl)+pc,r{k, uj;T = 0) (4.25) 



where ao^r{k,U!) = ZtTo(fc,^) is given by the expression H4.15|l but with the coupling constant replaced by the renor- 
malized coupling Zg^ = 

The continuum contribution is given by 



Pc,r{k,uj;T = 0) = - 



(To,r(Q^ 



(4.26) 



with 



1 f°° 



+M2) 



dP' 



(4.27) 



where we have made explicit the two particle threshold in the lower limit of the integral. 

The exact expression for Z given by the sum rule (|3.40() coincides with Z given by eqn. 1)4. 24(1 to lowest order in 
perturbation theory {(D{g'^)). 

Up to 0{g^) we can neglect cro(fc,Ci;) as well as TiQ^r{k,ijj) in the denominator of the continuum contribution (|4.26l) 
because > + M|) > mp and the denominator is never perturbatively small. Therefore to leading order in 
the coupling we can approximate 



Pc,r(«:,t^; J = 0) ~ - 



IT (g2-m|)2 • 



(4.28) 



The renormalized spectral function at finite temperature can be separated into the contributions from the different 
multiparticle cuts, 



Pr{k, uj, T) = /9/,r(fc, t^, T) + piiAk, ^1 T) 
where the contribution with support above the two particle cut is 



(4.29) 



[ao,rik,u!) + aa,rik,uj;T)] 



Pl,rik,UJ]T) = - 

[(g2-m|,)(l-no,.(Q2)j _ReS^,,(fc,c^;r)J + [ao,.(Q2) + a,^^(fc, c^; T) 
and that which has support below the light cone given by 



(4.30) 



pnAk,^\T) 



cTb{k,u;T) 



(Q2 _ ml)(l - UoAQ^)) - ReE^_^(fc,w;T) 



ab,rik,u);T) 



(4.31) 
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where again the renormalized quantities are obtained from the unrenormahzed ones by replacing g gr = Zg. 

Since pj r{k,uj) has support only for |a;| > i/fc^ + [Mf + M|) its denominator is never perturbatively small, there- 
fore to leading order 0{g'^) in the perturbative expansion it can be approximated by 

pUk, T) . i [ao.(fc o.) + a.,(fc, ^;T)] ^^ 3^^ 

TT (Cl!^ — — nip)'^ 

For pii^rik,uj) we must keep the full expression because for mp < (Mi — the denominator becomes pertur- 

batively small for uP' ^ k"^ + m?p. Therefore the final expression for the asymptotic distribution function H3.59|l to 
leading order in the coupling (0(g^)) is given by 

N{k,T)=^Ni{k;T)+NH{k,T) (4.33) 
where the different contributions reflect the different multiparticle cuts, namely 



k 

oc 



[1 + 2n{uj)] piiAk, LJ,T))dw -- (4.34) 



A/'/(fc,T) = - 



UJ 



n{uj) ao,r{k,uj) + aa,r{k,uj;T) + -aa,r{k,Lu;T)\ dcj , (4.35) 



2r!fc(c^2_f^2)2 

1 

where uJthik) = [k'^ + {AIi + Af2)^] ^ is the two particle cut. 

A. The approach to equilibrium: 

Before we study the asymptotic distribution function we address the approach to equilibrium. The time evolution 
of the (interpolating) number operator Nk{t) given by eqns. 1)3.1013.12(1 is completely determined by the real time 
evolution of the solution fk{t) given by eqn. H2.76|l . For mp < \Mi — M2I the particle acquires a width through 
the two body decay of the heavier particle in the bath and the particle pole is now embedded in the continuum for 
< (Ml — M2)^, which is the relevant part of the spectral density is (T;,(fc, cj, T) given in eqn. (|4.17() . In the Breit 
Wigner approximation, the spectral density is given by eqns. (|3. 1413. 1313. 15(1 with 

abik,Wk{T),T) 

r,(r)^z, ^^^^^^ (4.36) 

The real time evolution of the solution fk{t) in the Breit- Wigner approximation is given by eqn. ((3.20(1 . Figure 
Q) displays both the exact solution 1(2.76(1 and the Breit- Wigner approximation 1(3.20(1 for k = 0. The exact and 
approximate solutions are indistinguishable during the time scale of the numerical evolution as gleaned from this 
figure. 

The asymptotic long time evolution is determined by the behavior of the spectral density near the thresholds and 
is typically of the form of a power law '2^. However, such asymptotic behavior sets in at very long times, beyond the 
regime in which our numerical study is trustworthy. It is numerically exceedingly difhcult to extract the exponential 
relaxation from the power laws that dominate at asymptotically long time because the amplitude becomes very small 
in the weak coupling case. 

The main conclusion is that the distribution function approaches thermalization and becomes insensitive to the 
initial conditions for time scales t > = l/2Tk{T), where Tk{T) is the quasiparticle relaxation rate. 

B. The asymptotic distribution function: 

In the Breit- Wigner approximation and assuming a very narrow resonance near the physical particle pole 



I sign(cj) Ffc 



pn,r{k,uj,T) - - — ^ - sign(w)(5(w^ - nl) (4.37) 
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FIG. 7: The functions fk=o{t) and fk^{t) vs tmp for gV(167r^m|,) = 0.01 Mi = Amp ; M2 = mp ; T = lOmp. For these 
values of the parameters we find numerically: ZoiT) = 0.998 , Wo(r) = 0.973mp , ro(r) = 0.012mp. The exact solution and 
the Breit-Wigner approximation are indistinguishable. 



where in the second term on the right hand side the width has been neglected by assuming a very narrow resonance 
at Qk- Therefore in this narrow width approximation one would expect that the different contributions are given by 

J^ii{k,T)^n{nk) ; Afi{k,T)^0{g^) (4.38) 

where n{^k) is the Boltzmann distribution function for the stable particle. This rather simple analysis would lead to 
the conclusion that the corrections to the equilibrium abundance are perturbatively small. 

However, even for weakly coupled theories we expect this simple argument to be incorrect both in the high and 
low temperature regimes. The main reason for this expectation is that the approximation H4.37|l suggests that 
this argument neglects the fact that the spectral density has support for frequencies smaller than the position of the 
physical particle pole (namely for 7^ fifc). From the expression H4.34|l it is clear that the region of small w will lead to 
a substantial correction since for a; <C T the Bose-Einstein distribution function in 14.34|1 becomes n{Lu) ^ T/lo >> 1, 
thus the region of \uj\ < ilk and in particular \lu\ <^ T gives a non-trivial contribution to the abundance. The region 
of spectral density for > Qk will yield a much smaller, but non-negligible contribution. Furthermore in the high 
temperature limit T 3> fc,TOp,Mi.2 the width is expected to become large. This can be gleaned from the expression 
for ab{uj,k,T) in eqn. (|4.14|l . which for T ^ Wp'^ is proportional to T. This is clearly a statement that at high 
temperatures there is a large population of heavy particles which results in a larger number of processes Xi ^ X2 
in the medium, thereby increasing the width of the particle 4>. As the width of the spectral density near the physical 
particle pole increases, the spectral density has larger support in the small uj region, thereby increasing the off-shell 
contributions to the abundance. These arguments will be confirmed both analytically and numerically below. 

We now study numerically and analytically the asymptotic distribution function to assess precisely the magnitude 
and origin of the corrections to the equilibrium abundance. The parameter space is fairly large, thus we consider 
separately the cases of small momenta k ^ mp,Mi,M2,T and the case of large momenta k ^ mp, Mi, M2,T 
choosing the unit of energy to be the zero temperature pole mass of the particle, mp and keeping the value of the 
masses of the heavy fields fixed with Mi > M2 + mp. 

1. k = 

The limit fc = of the spectral density can be easily obtained from the expressions given above (|4.15I4.17|I . Of 
particular importance is the high temperature limit of (Tf,(0,w, T) since this contribution to the spectral density 
determines the width of the spectral function near the physical particle pole ro(T) given by eqn. (|4.3t)|) . 
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FIG. 8: The spectral density pii{k = 0,uj,T) vs uj/mp for g'/{16Tv'm%) = 0.01 Mi = 4mp ; M2 = mp and T/mp = 10,30, 100 
respectively. 



A straightforward calculation leads to the following result in the limit T ^ mp , Mi_2, 



ab{0,mp,T) g^T 


m% + (Af2 _ - 2ml{Ml + Mi) 


2 


2mp 87r2 


{Mf ~ MiY - m% 





(4.39) 



We note that this expression for the width is classical since restoring — > g^ft; T T/h the expression above 
is independent of h. This is a consequence of the fact that the high temperature limit is completely determined by 
the Rayleigh- Jeans part of the Bose-Einstein distribution function. As a result when the temperature is much larger 
than all mass scales, the width is proportional to T and the spectral density becomes wider, enhancing the off-shell 
contributions. Fig. JHJ displays the spectral density for several values of the temperature highlighting the broadening 
for large temperature. It is clear from this figure that at very high temperatures perturbation theory breaks down in 
this model since the width can become comparable to the physical mass or the position of the pole. This situation 
has been previously noticed in a scalar field theory at high temperatures, and a finite temperature renormalization 
group was introduced to provide a non-perturbative resummation[28j . 

Restricting ourselves to the regime in temperature within which perturbation theory is still reliable, namely for 
ro(T') <^ mp, we study the departure of the distribution function from the Bose-Einstein form (for fc = 0) numerically. 

Figure ^ displays the quantity 



A(T) ^ -^(fc^O, T) n(mp) ^^^^^^ 
n[mp) 

for a weakly coupled case in the range of temperatures 1 < T/mp < 20 for AIi — Amp , M2 = mp within which we 
find numerically that TQ(T)/mp < 0.1 which we use as a reasonable criterion for the validity of perturbation theory 
(see fig. IHl). 

This figure clearly indicates that even within the high temperature regime where perturbation theory is reliable and 
the spectral density still features a rather narrow Breit-Wigner peak, there are substantial departures from the Bose- 
Einstein form in the equilibrium distribution function. At low temperatures fig. @ clearly displays an exponential 
suppression and the distribution function essentially becomes the Bose-Einstein distribution. In this limit the width 
is extremely small and the spectral density is almost a delta function on the physical particle mass shell, and the 
off-shell effects are perturbatively small. 



2. A;>r,mp,Mi,2 



In the limit of large momenta several interesting features emerge: i) the width of the spectral density becomes 
very small, this is a consequence of the fact that there are very few heavy states for large momenta in the heat 
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FIG. 9: The ratio A(r) = (A/'(0, T) - n{mp))/n{mp) vs mp/T for / (WTr'^ml) = 0.01 Ah = 4mp ; Ah = mp. 



bath if the momentum is large. The width as a function of k is depicted in fig. H10() . which displays clearly this 
behavior, ii) As a function of the variable oj, the position of the peak in the spectral density becomes closer to the 
threshold for k ^ mp, Mi, 2. As a result of both these effects the spectral distribution becomes strongly peaked near 
threshold and the threshold moves to larger values of the frequency, thus leaving behind a larger region of the spectra 
off-shell for frequencies smaller than the position of the peak. The spectral density while small away from the peak is, 
however, non-vanishing and the fact that there is now a larger region in frequency lo below the (narrow) peak, brings 
about a competition of scales as can be understood from the following argument. The very narrow peak (almost a 
delta function a,t cu = ^Ik ^ k) leads to a contribution J\fjj{k,T) n{ilk), which for fc ^ T is <C 1. This on-shell 
contribution competes against the off-shell contributions from integrating the spectral density for to < flk which is 
also very small because ab{k,uj,T) /ftl <C 1 but for w ^ T is multiplied by the Bose enhancement factor ~ T/uj. The 
competition between the "on-shell" contribution n{Qk) and the off-shell contributions is studied numerically. 
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FIG. 10: The width r{k,T) in units of mp vs k/mp for fc > T for T = mp for /{16Tv'^m%) = 0.01 Ah = 4mp ; AI2 = mp. 



Fig. (|ll|l displays both the Bose-Einstein distribution function n(ilfc) and the asymptotic distribution function 
Af{k, T) in the limit fc > T, mp, Afi,2. 

It is clear from this figure that while the distribution function Mik, T) is strongly suppressed for fc 3> T it is larger 
than the Bose-Einstein distribution. The main reason for this enhancement is precisely the competition mentioned 
above, namely the position of the peak in the spectral density moves towards threshold which for large k corresponds 
to large values of the frequency lu. Therefore there is a large region in which the spectral density is very small but 
non- vanishing for lu < Q^- Clearly the part of the spectral density with support for uj > yields a much smaller 
contribution to the distribution function. Furthermore, for uj ^ T the factor n(uj) ^ T/lu 3> 1 which enhances further 
the off-shell contributions. 
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FIG. 11: The distribution functions M{k,T) given by eqn. 14.3311 (solid line) vs. the Bose-Einstein distribution NBE{k,T) 
(dashed line) as a function of k/mp for T = mp; = 0.01; Ml = 4mp ; M2 = mp. 



These results in the different regimes can be summarized as follows: 

• In the high temperature regime the larger abundance of heavier particles in the bath leads to a broadening of 
the spectral density. This broadening in turn results in a larger off-shell contribution to the abundance A/'(fc, T) 
and an enhancement of the distribution function over the Bose-Einstein result. The off-shell region of small 
frequency yields a substantial contribution because of the factor n{uj) ~ T/u in H3.59|l . In the model considered 
perturbation theory breaks down at high temperature and the imaginary part on-shell becomes classical. This 
situation is akin to the case of a self- interacting bosonic field theory studied in ref. [2^. A high temperature 
renormalization group resummation program such as in ref.|28| may be required to provide a non-perturbative 
resummation. 

• For momenta much larger than the mass scales and the temperature there is also a large enhancement of the 
distribution function JV{k, T) over the Bose-Einstein result. In this case the spectral density features a very 
narrow resonance near the position of the physical pole aX uj ^ U.^^ which however moves closer to threshold. 
For large k the off-shell region of support of the spectral density becomes larger and though the spectral density 
is strongly suppressed, the off-shell contribution from the region lo < competes with the contribution from 
the on-shell pole, namely the Bose-Einstein distribution function n(ilfe) because for k » T n{Q,k) <C 1. The 
off-shell contribution from the region uj << ^1^ is comparable to or larger than n(f2fe) for k ^ T and is enhanced 
in the region w <C T by the factor n(Lu) ^ T/lu. 

While these results may be particular to the model studied, we expect most of these features to be robust and 
fairly general. In particular at high temperature it is physically reasonable to expect a thermal broadening of the 
spectral density either from collisions, many-body decays or Landau damping as in the case studied here. Broadening 
of the spectral function yields a larger contribution from the small tu region which is enhanced further by the factor 
n{uj) ^ T/lo for w <C T. Therefore a substantial departure of the distribution function A/'(fc, T) from the Bose-Einstein 
distribution is expected at high temperature. A possible breakdown of perturbation theory in the high temperature 
regime may require the implementation of a non-perturbative resummation procedure akin to that introduced in 
ref.|28|. At low temperatures, much lower than the mass and momentum scales a departure from simple Bose- 
Einstein is also expected. In this case even though the spectral function features a sharp and narrow peak at a 
position very near the physical particle pole, the Bose-Einstein distribution function is very small. Hence the off-shell 
region oj << 17^ of the spectral function will lead to a substantial contribution which is further enhanced by the 
factor n{uj) ^ T/lo for oj <^ T. Again because the temperature is much smaller than any of the scales, the spectral 
density will be exponentially suppressed off shell and the equilibrium abundance will reflect this suppression, but just 
as in the case studied here, may still be larger than the simple Bose-Einstein abundance. Of course our study within 
this particular model serves only as a guidance and the details of the enhancement will depend on the theory under 
consideration, but the main lesson learned here is that the off-shell, small frequency region of the spectral density 
yields a substantial contribution to the equilibrium abundance in interacting theories. 
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V. BOLTZMANN KINETICS IN RENORMALIZED PERTURBATION THEORY 



It is important to understand the origin of the differences between the quantum kinetic equation for the distribution 
function H3.10|l and the usual quantum Boltzmann equation. Therefore in this section we provide a derivation of the 
quantum Boltzmann equation in renormalized perturbation theory to highlight the origin of the different equilibrium 
abundances. We assume that the bath is in equilibrium just as we did in our derivation of the effective action and 
the time evolution for the distribution function in the previous sections. 

The quantum Boltzmann equation is a differential equation for the single particle distribution function. However, 
as we have discussed in detail above, the physical particles have mass mp and the Heisenberg field operators create 
physical particles out of the vacuum with an amplitude determined by the wave function renormalization. Therefore 
in order to account for the mass and wave function renormalization, and to obtain the kinetic Boltzmann equation 
for the physical particles it is convenient to re-write the Lagrangian by introducing counterterms, namely 



E 



i=l 



9r^r Xl X2 + Ccount + A«t [Xl 5 Xs] (5.1) 

(5.2) 



where — vZg and we assume that the renormalization aspects of the fields xi.2 had already been included in 
Cintixi 1 X2] • The counterterms in 1)5.2(1 are treated systematically in perturbation theory along with the cubic 
interaction. Note that Z — I , Am^ are both of 0{g^). 

The renormalized field is expanded in terms of creation and annihilation operators of physical Fock states, 



k 

mh 



— k 



(5.3) 
(5.4) 



and a similar expansion for the bath fields xii X2 in terms of creation and annihilation operators and the corresponding 

(1 2) 

frequencies ■ The total interaction Lagrangian is 



Lint - E E Xi.pXa -p--fe + Lc 



(5.5) 



k P 



where Lcount is the counterterm Lagrangian. The kinetic Boltzmann equation for the occupation number of the Fock 
quanta of the field <i> is 



dt 



dNk 
dt 



dNk 
dt 



(5.6) 



gain 



loss 



The gain and loss terms are obtained by calculating the transition probabilities per unit time for processes that lead 
to the increase (gain) and decrease (loss) the occupation number, namely dNk{t)/dt — dP^/dt. Within the framework 
of the kinetic description such calculation is carried out by implementing Fermi's Golden rule. The processes that 
lead to the increase or decrease of the population are read-off the interaction and energy conservation emerges as 
a consequence of taking the long time limit as is manifest in Fermi's Golden rule. The cubic interaction term in 
Lint gives rise to several different processes which are gleaned by expanding the product in terms of the creation 
and annihilation operators of all the fields involved. The different phases that enter in such terms determine the 
energy conservation delta functions in Fermi's Golden rule. Some of the processes are depicted in fig. Q. When 
mp < Mi,M2 the quanta of the field $ cannot decay into those of the bath fields, however if Mi > M2 + mp (or 
M2 > Ml +mp) the heavier bath field can decay into particles of $ therefore increasing the population. This process 
is depicted in fig. ©(b)). The inverse process contributes to the loss term. Let us consider the case Mi > M2 + mp 
(the case M2 > Mi + mp is similar by M2 Mi). The only process that leads to the gain in the population by 
energy conservation is xi ^ ^ X2 and consequently the only process that leads to the loss of population with energy 
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conservation is the inverse process of annihilation $ X2 — > Xi- The calculation of the gain and loss terms is as follows: 
consider the initial Fock state |7Vg, n^^\ n*^^-) where N is the occupation of particles of $ and n^^'^^ that for the 
respective bath fields. To lowest order in the coupling g the interaction from the counterterm Lagrangian does not 
contribute to the gain or loss, but only to forward scattering since these terms are already of 0{g^) and the transition 
probabilities will be at least of (g^)- Thus to lowest order 0{g'^), the gain term arises from the following matrix 
element 



igr 



1 1 1 2 
1,12 -. — 1, 711. - 
' P ' p+k 



1 



p ' p+k' 



(5.7) 



The limits in the time integral had been extended to ±00 according to Fermi's Golden rule which leads to energy 
conservation. The calculation of this matrix element is straightforward, taking the absolute value squared of this 
matrix element, summing over p and averaging over the occupation numbers of the particles in the bath, which is 
assumed in equilibrium, one obtains the inclusive transition probability 



gain 



2Vli. 327 



p+k P 



(5.8) 



f 1 2) 

where n(u; -. ' ) are the Bose-Einstein distribution functions since the thermal bath is assumed to remain in equilibrium. 



27t5 (nk+ wl^^- - oj^V) = 2tt5 (Vlk+ oj^?^ - ^1 t where t is the 

\ p+k P J \ p+k P J 



To obtain the above expression we have used 

total interaction time. A similar calculation leads to the total transition probability for the loss process: 



P - 2* AT 

^loss — 7:7^ or. O ^^k 



(1) (2) 
P p+k 



2nk 327r2 

The kinetic equation can now be written in the following form 



l + n(4^))] n{oj^l.)s{nk + 



2 1 
UJ - — uil 
p+k P 



(5.9) 



dNk 
dt 



(5.10) 



The gain and loss rates are given by 
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P p+k 



d^p 



+ n(.«)] niJ^^si^k+u:^,-^') 



(5.11) 
(5.12) 



Since the bath particles are in thermal equilibrium with a Bose-Einstein distribution function the detailed balance 
condition follows, namely 



(5.13) 



The solution of the Boltzmann kinetic equation H5.10|l is the following 



Nk{t)^n{^lk) + [Nk{'^)-n{Slk)]e 



-7fc * 



(5.14) 



where 



7fc = Ik - 7fc = 



2Vtu 327r2 



d^p 



1 (2) 
P p+k 



/ (2) s / (1) J r (f^ (1) , (2) 



(5.15) 
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Comparing this expression with those for the imaginary part of the self energy given by 114.1114.1214.1314.1^ it is 
straightforward to see that 

where 

ImSf (fc, fifc, T) = <7bAk, ^k, T) (5.17) 

This expression for the relaxation rate should be compared to the decay rate for the single quasiparticle Tk given 
by eqn. (|3.16I3.20|I ^. Since the quasiparticle residue in perturbation theory is Zk{T) = 1 + 0{g^) and the difference 
between the quasiparticle frequency yVfc(T) and the single particle frequency flk is of 0{g'^) to leading order in the 
coupling g, the relaxation rate of the distribution function jk and that of the single quasiparticle Tk (see eqn. ()3.16f) ') 



7fc = 2rfc+0(/) (5.18) 

We have provided this derivation of the usual quantum Boltzmann equation and its solution in the case when the 
bath remains in equilibrium to highlight the similarities and differences with the real time evolution of the distribution 
function given by eqn. H3.10(l : 

• The derivation above clearly shows that the Fock states that enter in the matrix elements H5.7(l are the asymptotic 
free field Fock states associated with physical particles of mass mp. This is similar to the definition of the 
interpolating number operator (|3.2|l which is based on the free field asymptotic physical states, and includes 
both mass and wave-function renormalization. 

• By implementing Fermi's golden rule, namely taking the time interval to infinity, thereby enforcing the on 
shell delta function, extracting the linear time dependence and dividing by time to provide a local differential 
equation for the time evolution of the distribution function all memory aspects have been neglected. Namely 
implementing Fermi's golden rule results in neglecting memory effects, which in turn results in only on-shell 
processes contributing to the rate equation. Contrary to this, the real time evolution of the distribution function 
(|3.1U|I includes memory effects as is manifest in the time integrals H3.11I3.12|1 in (|3.1U|) . In turn these time 
integrals keep memory of the past time evolution, and at asymptotically long time lead to the full spectral 
density as manifest in eqn. (|3.23|) . not just an on-shell delta function. The presence of the full spectral 
density in the asymptotic distribution includes the off-shell contributions discussed in the previous section. This 
discussion brings to the fore that one of the main origins of the differences can be traced to memory effects and 
the fact that the real time evolution of the distribution function H3.1U|I is non-Markovian. The memory of the 
past time evolution translates in off-shell processes through the full spectral density. 

• As emphasized in section (|lHl the expression (I3.10|) for the quantum kinetic distribution function implies a Dyson- 
like resummation of the perturbative expansion and includes consistently the renormalization aspects associated 
with asymptotic single particle states, namely the correct pole mass and the wave function renormalization. 
The dependence of the asymptotic distribution function on the full spectral density is a consequence of the 
fluctuation-dissipation relation. 

VI. CONCLUSIONS AND DISCUSSION 

Motivated by a critical reassessment of the applicability of Boltzmann kinetics in the early Universe, in this article 
we studied the abundance of physical quanta of a field $ in a thermal plasma by introducing a quantum kinetic 
description based on the non-equilibrium effective action for this field. We focused on understanding the equilibrium 
abundance of particles that are stable in the vacuum and interact with heavier particles which constitute a thermal 
bath. 

The non-equilibrium effective action is obtained by integrating out the heavy particles to lowest order in the coupling 
of the field <i> to the bath but in principle to all orders in the coupling of the heavy fields amongst them. We show that 
the non-equilibrium effective action leads to a Langevin stochastic description with a Gaussian but colored noise and 
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a non-Markovian self-energy kernel. The correlation function of the noise and the non-Markovian self-energy kernel 
are related by a generalized fluctuation dissipation relation. The correlation functions are determined by the solution 
of this Langevin equation which furnishes a Dyson resummation of the perturbative expansion. We introduced a 
definition of the single physical particle distribution function in terms of a fully renormalized interpolating Heiscnberg 
number operator based on asymptotic theory. The real time evolution of this single particle distribution function is 
completely determined by the solution of the Langevin equation. 

We show that in a heat bath at finite temperature this number operator becomes insensitive to the initial conditions 
after a time scale ~ l/2rfc(r), where Tk{T) is the single quasiparticle relaxation rate. We prove that the asymptotic 
long time limit of this distribution function describes full thermalization of the $ particle with the thermal bath. The 
equilibrium distribution function depends on the full spectral density and includes off-shell corrections as a result of 
the non-Markovian real time evolution (with memory) and the fluctuation-dissipation relation. Its expression is given 
by eqn. 1)3. 59|) . We argue that while we obtained the distribution function in the case of a field linearly coupled to 
a thermal bath of heavier particles, the final form of the distribution function at asymptotically long time is much 
more generally applicable. 

In order to provide a detailed assessment of novel specific features of the distribution function in particular departure 
from the usual Bose-Einstein distribution, we considered a model in which the thermal bath is described by two heavy 
bosonic fields xi,2 coupled to the field $ as g xi X2, with Mi > AI2 + rrip and the pole mass of the field We 
obtained the real time effective action at one loop level. We find that the in-medium processes of two body decay of 
the heavier particle, and its recombination, namely xi ^ X2^ results in a width for the <I>-particle and a broadening 
of its spectral density. A detailed study of the single (physical) particle distribution function reveals substantial 
corrections to the Bose-Einstein distribution at high temperature as well as low temperature but large momentum. At 
high temperature the spectral density broadens dramatically and the off-shell contributions become very substantial 
resulting in an enhancement of the abundance with respect to the Bose-Einstein distribution. We found that at 
very high temperatures, perturbation theory breaks down and a resummation of the perturbative expansion via the 
renormalization group at finite temperature may be requiredH. This case must be studied further. 

In the limit where the momentum of the particle is much larger than the temperature and the masses, our analysis 
also reveals a substantial departure from the Bose-Einstein distribution. In this case the spectral density is sharply 
peaked near the (zero temperature) physical pole mass, but the position of the peak moves to higher energies. As a 
result, the spectral density features off-shell contributions in a large region of frequencies smaller than the position 
of the peak. The small frequency region is further enhanced by temperature factors and these off shell contributions, 
while exponentially small, compete with the exponentially small on-shell contribution which yields the Bose-Einstein 
distribution. As a result the distribution function, while strongly suppressed at high momenta much larger than 
the temperature (and mass scales), is considerably larger than the Bose-Einstein abundance predicted by the usual 
Boltzmann equation. 

In order to highlight the origin of the enhancement, we derived the Boltzmann equation in renormalized perturbation 
theory up to the same order in the coupling to the bath as the non-equilibrium effective action, which is the basis for 
the quantum kinetic description. This derivation makes manifest the origin of the discrepancy: the usual Boltzmann 
equation is based on Fermi's golden rule, which requires taking a long time limit in the transition probability. In 
taking the long time limit and extracting the asymptotic behavior of the transition probability energy conservation 
is manifest as an on-shell delta function, and all memory effects have been neglected. Furthermore in considering 
the transition probability in a gain-loss balance equation, interference phenomena have been neglected. As a result 
the Boltzmann equation neglects off-shell contributions. Precisely these off-shell contributions from the support of 
the spectral density away from its peak and near the particle mass shell, are responsible for the departure from the 
Bose-Einstein result. The enhancement over the Bose-Einstein distribution is a consequence of the off-shell support 
of the spectral density at frequencies smaller than the position of the peak. 

Although these results were obtained within the particular specific model studied here, the origin of the discrepancies 
suggests these to be much more general. The spectral density of a particle that is stable at zero temperature features 
an on-shell delta function below the multiparticle thresholds. However in a medium this peak will be broadened by 
different processes and the particle becomes a quasiparticle. This unavoidable feature of an interacting particle in a 
medium results in a broader spectral density with a region of support at frequencies smaller than the position of the 
peak, which leads to a larger contribution to the abundance as compared to the Bose-Einstein distribution which is 
the "on-shell" result. 

Cosmological consequences: An important feature of the distribution function (|3.59|) is that it is exponentially 
suppressed at low temperatures since all the intermediate states are heavy and therefore exponentially suppressed at 
low temperatures. Therefore the off-shell contributions are strongly suppressed leading to the conclusion that the low 
temperature abundance is exponentially suppressed. This is in agreement with the results of refs.0, Q. Therefore 
we do not expect the low temperature enhancement of the abundance to be of any practical relevance for cold dark 
matter candidates. 
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The consequences for the cosmic microwave background depend on the temperature regime. For temperatures 
much larger than the electron mass the photons in the plasma propagate as in-medium quasiparticles of two species: 
longitudinal and transverse plasma excitations (plasmons). The plasma frequency in the high temperature regime 
is of the order ujpi oc ^/a^r„ T|23l |2^ . The corrections to the dispersion relations (plasma frequency) arise from 
intermediate states of electron-positron pairs and yield a contribution to the spectral density with support below 
the light cone. These are Landau damping processesjl^ while those that yield the width arise from Compton 
scattering and pair annihilation and are of higher order. The plasmon width (up to logarithmic corrections) is of 
order F cx a^^T. Thus the spectral function for pho tons features support both above and below the light cone, the 
latter is a result of Landau damping processes |2a l23 | . This latter contribution is important because it yields support 
in the small frequency region which is Bose enhanced. Both the plasma frequency and the width are strong functions 
of temperature and we expect substantial corrections to the power spectrum of the cosmic microwave background for 
T ^ 1 Mev. However, these potential corrections are observable only indirectly, possibly through nucleosynthesis. 
For temperatures well below the electron mass the lowest order 0{aem) correction to the spectral density arises from 
an electron loop and an electron-positron loop (we ignore the contribution from protons). The former gives a Landau 
damping cut below the light cone akin to the contribution (14.1411 ^2^, and the latter gives a two particle cut above 
the pair-production threshold. Both are off-shell contributions and yield corrections to the spectral density which are 
proportional to the electron number density ( equal to the proton number density) rig ^ Xe{^bhQ){l + z)'^ x lO^^cm"'^, 
with Xe the ionization fraction. The width of the spectral density near the mass shell results from Compton scattering 

and is of order Ofg^. It is approximately given by F ~ aTne\J~^^ and ar is the Thompson scattering cross section. 
During recombination the ionization fraction diminishes precipitously within a window of redshift Az ~ 100 which is 
the width of the last scattering surface. This rapid vanishing of the ionization fraction and consequently of the (free) 
electron density entails that the broadening of the spectrum and the spectral distortions become vanishingly small at 
the end of recombination. At decoupling the mean free path is comparable to the size of the horizon and the spectral 
density for photons is basically that of free field theory. Hence recombination erases any observable vestige of spectral 
distortion through many body processes and spectral broadening, thus there are no observable consequences of these 
effects in the CMB. 

However we expect that our results may be potentially relevant in the high temperature limit for the kinetics of 
baryogenesis in the Early Universe. We expect to address these issues in future work. 
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APPENDIX A: CALCULATION OF THE IMAGINARY PART OF THE SELF-ENERGY 

The imaginary part of the self-energy is given in the text, eqn. H4.9|) . 

Integrating over dpo, dqo and then performing the transformation p —p— k in all the integrals involving n{uj\ -), 
we can write 

ImS^(w, k)^ao + ai + an + {cj^i}i ~ cjf},) + (frj;) - af^) (Al) 

where 
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Obviously, cto represents the zero temperature contribution. Note that an and ct^^ can be obtained by exchanging 
Ml and M2 in aj and cr^jj respectively. Thus, we will only outline the main steps in computing do, cri and crj^j in 
this appendix. First of all, let tOp = ojp^^ and z = i^-^j:- Then, we have 

ctq + cr/ = — — sign(a;) / [\ + n{ujp)]dujp I 5{\uj\ - ujp - z) dz (AT) 
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where 



Ml 



'2 ± 2fc^C^2 _ ^J2 ^ ^2 _ (^^2 _ ^,^2)^ 

Without loss of generality we can assume that Mi > M2 for convenience. For the integral to be non- vanishing, we 
require that 

z^ < z = \u}\ - ojp < z+. (AlO) 
Squaring both sides twice properly, these two inequalities can be reduced to fitOp) < where 

f{cjp) = 4{\Ljf ~ fc2)^2 _ 4|^|(|^|2 _ ^ (|^|2 _ ^)2 ^ 4^^^2 

and a = k'^ ~ {Ml — Notice that the graph f{u!p) against ujp represents a conic with positive y-intercept. 

Solving f{LOp) — for LUp, we obtain 



± _ |w|(|w|2 - a) ± fcv/(|t^P-a)2-4(|c^|2-fc2)M2 



"^"< = 2(|..|2-P) (^12) 

There are two possibilities: (i) jcjp — /c2 ^ (^y-j ^2_ |^|2 ^ q p^j. ^2_ |^|2 ^ graphs with /(wp) against ujp show 
that condition (jAlOp can be satisfied only if ujp > uj^ but algebraic calculation indicates that |u;| — w~ < 0. Thus, 
condition (| A 10(1 can never be satisfied and this solution should be ignored. For jwp — k"^ > 0, we have \uj\'^ — a > 0. 

A detailed analysis of /(wp) as well as and \uj\ — ujp as functions of uOp reveals that that condition IjAlOp can 
always be satisfied for LOp < ujp < uj^ and > ■\/fc2 + M| + Mi. For the discriminant in ui^ to be positive, we 
require that \uj\ > y^k^TWi^^W or \lu\ < ^k^ + {Mi- M2Y. Since ^/WTM^ + Mi > y/k^ + (Mi - 1/2)2, 
we can only pick up \uj\ > ^Jk"^ + [Mi + 1/2)2. As a result, we conclude that 
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^0 = T|^signHe[|c^|2-P-(Mi + M2)2](c^+-c^;), (A13) 



= ^ sign(.) e[ - e - (M, + M,f ] In i^^^-^] • (AM) 



Now, we proceed to compute crf^f. 



n{ujp) dujp / S{uj + ujp — z) dz. (A15) 

Ml 



For the integral to be non- vanishing, we require that 

z~ < z — oj + LUp < z'^ {A16) 

which can be reduced to g{LOp) < where 

g{LUp) = 4(w^ - + 4a;(tj^ - a)wp + (w^ - a)^ + Akruf. (A17) 

Solving g^LOp) = for tjp, we obtain 



= i^) = 2(c.2 - P) ^^^^^ 

First, note that — > ± A; as — *■ oo. Then, drawing graphs with g{uOp) against ujp and diagrams with 
and w + Wp against cjp, we observe that condition ljA16|l is always satisfied for k"^ — lo^ > with Wp > ^~{ui). For 
uj'^ ~ k'^ > 0, we have |wp — a > and graphs with g{cOp) against uip show that condition ljA16|l can be satisfied only 
if ti) < and < Wp < Moreover, an algebraic calculation indicates that both u) + ^~ < and uj + S,^ < unless 
Lu^ — k'^ < Ml — A/| . Additionally, for the discriminant in to be positive, we require that oJ^ — k^ > {Mi + 
or _ fc2 ^ (-^^ _ Afj)^. The condition uo"^ ~ k'^ > (Mi + Af2)^ contradicts - fc^ < M^^ - M|. Hence, we must 
take — fc^ < (Ml — M2)^. Graphs of and w + Wp against ojp confirm that condition (|A16p is always satisfied for 
< - A:^ < (Ml - M2f. As a result, we have 

a« - afA = ^ e(fc^ - o.^) In ( I^Jl!^\ (A19) 

^\£kp sign(^) ©('^^ - k^) 0[fc' + (Ml - M2f - ] In (A20) 

(A21) 

where Wp are the roots given by (|A12I) . For k^ — u"^ > and uj > 0, ^~{—uj) = |a;~| and S,p{uj) — \oJp\- For 
/c^ — > and w < 0, ^~{—lu) = and ^p'{oj) = \iLJp\. Therefore, we conclude that 



^1// - <^?ii = Bign(^) e[ fc^ + (Ml - M^r - ] In (l-^^^ . (A22) 

Finally, to obtain cr^y — afy , we simply need to exchange Mi and M2 in — afjj . 
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